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Abstract 



In the present work we discuss inflation, dark matter and cosmological evolu- 
tion in the context of the brane-world scenario. Being string theory inspired, 
the brane-world models provide corrections to the General Relativity, which 
is considered to be the low-energy limit of string theory. We find that novel 
cosmologies are obtained, which potentially can provide answers to some 
of the longstanding problems of modern cosmology, such as the origin and 
nature of dark energy. At the same time the successes of the standard four- 
dimensional cosmology are preserved, and in some cases the treatment in the 
framework of brane cosmology is even more satisfactory. 
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Chapter 1 



Introduction 



The standard electroweak model is a mathematically-consistent renormaliz- 
able field theory which predicts or is consistent with all experimental facts. 
It successfully predicted the existence and form of the weak neutral current, 
the existence and masses of the W and Z bosons, and the charm quark, 
as necessitated by the GIM mechanism. The charged current weak interac- 
tions, as described by the generalized Fermi theory, were successfully incor- 
porated, as was quantum electrodynamics. When combined with quantum 
chromodynamics for the strong interactions and general relativity for classi- 
cal gravity, the standard model is almost certainly the approximately correct 
description of nature down to at least 10~^^cm, with the possible exception 
of the Higgs sector. However, the theory has far too much arbitrariness to 
be the final story. For example, the minimal version of the model has 21 free 
parameters, assuming massless neutrinos and not counting electric charge 
(Y) assignments. Most physicists believe that this is just too much for the 
fundamental theory. The complications of the standard model can also be 
described in terms of a number of problems. 



1.1 Gauge Problem 

The standard model is a complicated direct product of three sub-groups, 
SU{3) X SU{2) X f/(l), with separate gauge couplings. There is no explana- 
tion for why only the electroweak part is chiral (parity- violating). Similarly, 
the standard model incorporates but does not explain another fundamental 
fact of nature: charge quantization, namely why all particles have charges 
which are multiples of e/3. This is important because it allows the electrical 
neutrality of atoms {\qp\ = \qe\)- Possible explanations include: grand uni- 
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tied theories, the existence of magnetic nionopoles, and constraints from the 
absence or cancellation of anomalies. 



1.2 Fermion Problem 

All matter under ordinary terrestrial conditions can be constructed out of the 
fermions (z/g, e~, -u, d) of the first family. Yet we know from laboratory studies 
that there are > 3 families: {u^j,, fi~ ,c,s) and {urjT^ ,t,b) are heavier copies 
of the first family with no obvious role in nature. The standard model gives 
no explanation for the existence of these heavier families and no prediction 
for their number. Furthermore, there is no explanation or prediction of the 
fermion masses, which vary over at least five orders of magnitude, or of the 
CKM mixings. There are many possible suggestions of new physics that 
might shed light on this, including composite fermions; family symmetries; 
radiative hierarchies, in which the fermion masses are generated at the loop- 
level, with the lighter families requiring more loops; and the topology of 
extra space-time dimensions, such as in superstring models. Despite all of 
these ideas there is no compelling model and none of these yields detailed 
predictions. The problem is just too complicated. Simple grand unified 
theories don't help very much with this, except for the prediction of mi, in 
terms of m,- in the simplest versions. 



1.3 Higgs/hierarchy Problem 

In the standard model one introduces an elementary Higgs field into the 
theory to generate masses for the W, Z, and fermions. For the model to be 
consistent the Higgs mass should not be too different from the W mass, Mfj = 
0{M^). If Mh were to be larger than Myy by many orders of magnitude 
there would be a hierarchy problem, and the Higgs self-interactions would be 
excessively strong. Combining theoretical arguments with laboratory limits 
one obtains Mh ^ 1 TeV. 

However, there is a complication. The tree-level (bare) Higgs mass re- 
ceives quadratically-divergent corrections from loop diagrams. One finds 

Mf, = (M^)barc + 0(A, /, h^)A^ (1.1) 

where A is the next higher scale in the theory. If there were no higher scale one 
would simply interpret A as an ultraviolet cutoff and take the view that Mh 
is a measured parameter and that (M//)barc is not an observable. However, 
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the theory is presumably embedded in some larger theory that cuts off the 
integral at the finite scale of the new physics. For example, if the next scale 
is gravity then A is the Planck scale Mp = G^ ~ 10^^ GeV. If there is a 
simple grand unified theory, one would expect A to be of order the unification 
scale Mx ~ 10^^ GeV. Hence, the natural scale for Mh is 0(A), which is much 
larger than the expected value. There must be a fine-tuned and apparently 
highly contrived cancelation between the bare value and the correction, to 
more than 30 decimal places in the case of gravity. If the cutoff is provided by 
a grand unified theory there is a separate hierarchy problem at tree-level. The 
tree-level couplings between the Higgs field and the superheavy fields lead to 
the expectation that Mh is equal to the unification scale unless unnatural 
fine-tunings are employed. 

One solution to this Higgs/hierarchy problem is the possibility that the 
W and Z bosons are composite. However, in this case one would apparently 
eliminate the successes of the SU{2) x f/(l) gauge theory. Another approach 
is to eliminate elementary Higgs fields in favor of a dynamical mechanism 
in which they are replaced by bound states of fermions. Technicolor and 
composite Higgs models are in this category. The third possibility is super- 
symmetry, which prevents large renormalizations by enforcing cancellations 
between the various diagram contributions. However, most grand unified 
versions do not explain why (Mw/Mx)'^ is so small in the first place. 



1.4 Strong CP Problem 

Another fine-tuning problem is the strong CP problem. One can add an 
additional term ^^glFF to the QCD lagrangian which breaks P, T and 
CP symmetry. F^y = e^jyQ/3F"'^/2 is the dual field. This term, if present, 
would induce an electric dipole moment d]x for the neutron. The rather 
stringent limits on the dipole moment lead to the upper bound O < 10~^°. 
The question is, therefore, why is so small? It is not sufficient to just 
say that it is zero because CP violation in the weak interactions leads to 
a radiative correction or renormalization of O by O{10~^). Therefore, an 
apparently contrived fine-tuning is needed to cancel this correction against 
the bare value. Solutions include the possibility that CP violation is not 
induced directly by phases in the Yukawa couplings, as is usually assumed in 
the standard model, but is somehow violated spontaneously. then would 
be a calculable parameter induced at loop level, and it is possible to make 
sufficiently small. However, such models lead to difficult phenomenological 
and cosmological problems. Alternately, becomes unobservable if there is 
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a massless u quark. However, most phenomenological estimates are not con- 
sistent with rriu = 0. Another possibihty is the Peccei-Quinn mechanism, in 
which an extra global f/(l) symmetry is imposed on the theory in such a way 
that O becomes a dynamical variable which is zero at the minimum of the 
potential. Such models imply the existence of very light pseudoscalar par- 
ticles called axions. Laboratory, astrophysical, and cosmological constraints 
allow only the range 10* — 10^^ GeV for the scale at which the above ^(1) 
symmetry is broken. 



1.5 Graviton Problem 

Gravity is not fundamentally unified with the other interactions in the stan- 
dard model, although it is possible to graft on classical general relativity by 
hand. However, this is not a quantum theory, and there is no obvious way to 
generate one within the standard model context. In addition to the fact that 
gravity is not unified and not quantized there is another difficulty, namely 
the cosmological constant. The cosmological constant can be thought of as 
energy of the vacuum. The energy density induced by spontaneous symmetry 
breaking is some 50 orders of magnitude larger than the observational upper 
limit. This implies the necessity of severe fine-tuning between the generated 
and bare pieces, which do not have any a priori reason to be related. Possible 
solutions include Kaluza-Klein and supergravity theories. These unify grav- 
ity but do not solve the problem of quantum gravity or yield renormalizable 
theories of quantum gravity, nor do they provide any obvious solution to the 
cosmological constant problem. Superstring theories unify gravity and may 
yield finite theories of quantum gravity and all the other interactions. It is 
not clear whether or not they solve the cosmological constant problem. 



1.6 Dark sector problem 

On the other hand Cosmology, the science of the Universe, has its own Stan- 
dard Model. It is the so-called Hot Big-Bang model. The expansion of the 
Universe, the Bing-Bang Nucleosynthesis and the Cosmic Microwave Back- 
ground Radiation have established the theoretical framework of Cosmology 
(which is based essentially on the Theory of General Relativity and the Cos- 
mological Principle, namely that the Universe is homogeneous and isotropic 
on large scales) into the Standard Model of modern cosmology. In the last 
decade or so we have entered an era of precision cosmology. The basic quan- 
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titles have been measured and now what remains Is to try to understand 
them. What Is surprising, Is the fact that the Universe today seems to ex- 
pand with an accelerating rate, while one would expect that because of the 
attractive nature of gravity the universe should be decelerating. This means 
that the dominant component In the universe Is some strange "material" 
with negative pressure, called Dark Energy. The rest of the universe consists 
of photons, neutrinos (only a tiny fraction of the energy budget) and non- 
relatlvlstlc matter, most of which Is In some unknown non-baryonlc form, 
called Dark Matter. So we see that In Its majority the Universe consists of 
something that we do not know what It Is. The nature and detailed charac- 
teristics of dark matter and dark energy are the major theoretical challenges 
for modern cosmology. 



1.7 The brane- world idea 

Recently It has been suggested that there might exist some extra spatial 
dimensions. Of course this Is not a new proposal. Instead It Is essentially a 
revival of the old Kaluza-Kleln Idea. In the traditional Kaluza-Kleln sense 
the extra dimensions are compactlfied on a small enough radius to evade 
detection In the form of Kaluza-Kleln modes. However now there Is a different 
setting where the extra dimensions could be large, under the assumption 
that ordinary matter Is confined onto a three-dimensional subspace, called 
brane (more precisely "3-brane", referring to the three spatial dimensions) 
embedded in a higher dimensional space, called bulk. In fact the idea that 
we might be living inside a defect embedded in a higher dimensional space 
has already a long history. In the context of an ordinary higher dimensional 
gauge field theory it was proposed that we might live on a codlmenslon one 
solltonlc object. However, it was soon realized that in contrast to scalar and 
spln-1/2 fields, it would be difficult to confine gauge fields on such an object. 

The situation is drastically different in the context of type-I string theory. 
A few developments led to an exciting possibility and renewed Interest in 
the whole idea. First, with the discovery of the D-branes as an essential 
part of the spectrum in type-I string theory, one could conjecture that we 
inhabit such a D-brane embedded in a ten-dimensional bulk. The usual 
solltonlc defect of field theory was thus replaced by an appropriate collection 
of D-branes, which by construction confine the gauge fields, together with 
all the ingredients of the standard model. All known matter and forces lie 
on our brane-world, with the exception of gravity which acts in the bulk 
space as well. It was however pointed out that the gravitational force on the 
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brane was consistent with all laboratory and astrophysical data as long as 
the extra dimensions were smaller than a characteristic scale. This led to 
the exciting possibility of two extra dimensions in the sub-millimeter range. 
Furthermore, it was demonstrated in the context of an appropriate effective 
five-dimensional theory of gravity, that once we take into account the back 
reaction of the brane energy-momentum onto the geometry of spacetime, the 
graviton is effectively confined on the brane and Newton's law is reproduced 
to an excellent accuracy at large distances, even with a non-compact extra 
dimension. 

The present work comprises my research in the field of brane cosmology. 
Brane-worlds open-up new ways to attack the Dark Matter, Dark Energy 
problems and these were exploited in what follows. Our work is organized as 
follows: In Part I, which includes chapters 2 and 3, we present the Theoretical 
framework of our discussion. We introduce, briefly and for completeness, the 
fundamentals of the Brane- World scenario, the Standard Models of particle 
physics and of Cosmology, as well as some basic knowledge about Supersym- 
metry. In Part II, the rest of the present work, we present our research in 
the field of brane cosmology. In chapter 4 we discuss axino dark matter in 
the brane-world scenario. We present sneutrino inflation and supersymmet- 
ric hybrid inflation in chapters 5 and 6 respectively. The role of brane-bulk 
energy exchange and a concrete model with our universe as a global late-time 
attractor is discussed in the seventh chapter. We summarize our results and 
finish with some conclusions in the last chapter. 



PART I 



The Theoretical Framework 



Chapter 2 



The brane- world scenario 



In the present section we shall describe various realizations of the brane- world 
idea [T]. 



2.1 Randall-Sundrum localization 

We will consider the case of a five-dimensional bulk with coordinates x^^ = 
(x^, y), fi = 0,1, 2, 3. We also consider a three-brane located at y = 0. Apart 
from the five-dimensional Einstein term we also have a constant energy den- 
sity (brane tension) on the brane, and a non-vanishing cosmological constant 
in the bulk. We can summarize the effective action as 

S = J dy d^x^{M^R - A) - / d^^dyS{y)^Vb (2.1) 

where gab is the induced metric on the brane gab = gMN^r%r- We will 
pick a static gauge for the brane coordinates C,"" = x"". To simplify matters, 
we will also consider an orbifold structure under y -^ —y. Thus, the two 
parts of space-time separated by the brane are mirror symmetric around the 
position of the brane. 

We would like to solve the equations of motion stemming from action 
(12.11) . Let us first seek solutions invariant under the orbifold action which 
are fiat along the brane, and depend only on the fifth coordinate y. 

The ansatz for the five-dimensional metric is 

ds^ = e^^^y^dx^'dx^ + dy^ (2.2) 

In order for the fiat-brane ansatz to provide a solution, the two vacuum 
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energies must be related 

V^ = -12AM3 (2.3) 

This implies that the vaccum energy A must be negative. We will also define 
the RS (AdS) energy scale 

^ = -^ (2.4) 

The gravitational equations have the solution 

A{y) = -K\y\ (2.5) 

The space on the one side of the brane is a slice of AdS^ patched-up with 
its mirror image at y = 0. Indeed defining r = e^"^ for y > and scaling 
x^ -^ x^ / K we obtain 

ds^ = -—^(dr^ + dx^dx^), r > 1 (2.6) 

which is the r > 1 slice of AdSr, in Poincare coordinates with AdS energy 
scale K. Note that the orbifold has removed the boundary of AdSr, at r = 0. 

An interesting further question concerns the effective interactions medi- 
ated by gravity in this background. To find them we must study the small 
fiuctuations around this solution. Direct variation of the equations along the 
brane longitudinal directions and gauge fixing gives a scalar equation for the 
static graviton propagator 

M%-e-'^Vl - dl - Aj^dy)G{x- y) = 6{y)6'^'\x) (2.7) 

This can be Fourier transformed along the a;* coordinates obtaining 

M'{-e-'^f - dl - AA'dy)G{f>- y) = 5{y) (2.8) 

Imposing the symmetry G{p,y) = G{p, —y) we obtain the solution 

G{y,p) = Bw'K,C^) (2.9) 

where p = |p|, w = e^'^'^', K2 is the standard Bessel function and 

We can now investigate the force mediated by the graviton fiuctuations 
on the brane by evaluating 

GiP^^) = 7T^ ^-^M{ (2-11) 



1 


K2 


(f) 


2M3p 


Ki 


(f) 
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The static gravitational potential between two unit sources on the brane 
(upon transforming back to configuration space) becomes 

V{r) = ^ r dppszn{pr)G{p, 0) = ^-^ + SV{r) (2.12) 

with 

where here r^ = x^ is the spatial distance on the brane. We can now compute 
the correction to the gravitational potential for two extreme cases. 

For Kr ^ 1 the main contribution to the integral ()2.13p comes from 
small q and we obtain 

Thus, at long distances gravity is four- dimensional with sub-leading correc- 
tions. The effective four-dimensional Planck scale is M^ = M^/K. 

For Kr <^ 1 the main contribution to the integral comes from large q 
and we obtain 

•^^W - , oL o (2-15) 

Thus, at short distances gravity is five-dimensional. This is completely anal- 
ogous to compactification with radius 1/K. The RS setup is thus an alterna- 
tive mechanism to compactification for transforming five- dimensional gravity 
into four- dimensional at large distances. 



2.2 Brane-Induced gravity 

We now describe an alternative realization of four- dimensional gravity that 
comes under the name of brane-induced gravity. First we consider the sim- 
plest case of a five-dimensional bulk with coordinates x^ = {x'^,y),fi = 
0,1,2,3 . We also consider a three brane located at y = 0. Apart from 
the five-dimensional Einstein term we would like to study the effects of a 
four- dimensional Einstein term localized on the brane. The relevant action 
is 

M^ / dyd^x^R + r^ I d^^dy6{y)^R (2.16) 

where Qab is the induced metric on the brane and R the induced Ricci scalar. 
We also parameterized the coefficient of the four- dimensional term in terms 
of a new length scale Vc- 

10 
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We are interested in the gravitational interaction, generated by the action 
()2.16|) . as perceived on the brane. We will evaluate first the static propagator 
of ()2.16|) . Although there is interesting physics in the tensor structure, we 
will neglect it for the moment and consider instead the scalar propagator. 
Placing the source on the brane (at the origin) we must solve 

M\yl + dl + r,5{y)Vl)G{x, y) = -6{y)6'^'\x) (2.17) 

Fourier transforming in the 3-spatial coordinates x we obtain 

M'{f-dl+r,6iy)f)Gip,y) = 5iy) (2.18) 

The solution can be found by first solving the equations away from the po- 
sition of the brane, and then matching along y = 0. The result is 

e-b1lyl 

For the source and the probe being on the brane y = the static propagator 
becomes 

G(p, 0) = ,,,, , ? ^ (2.20) 

^■^' ' M3(2|p1+rcp2) ^ ' 

By Fourier transforming back we obtain the static gravitational potential 

1 r°° 

V(r) = ^— / dppsin(pr)G(p,0) (2.21) 

We are now ready to study the behavior of the gravitational force in various 
regimes. 

For long distances pvc <^ 1 the potential can be approximated as 

^W - M^ (2.22) 

This is the behavior of five-dimensional gravity with Planck scale M. 
For short distances pvc ^ 1 we obtain 

V{r) ~ — ^ (2.23) 

This is the behavior of four-dimensional gravity with effective Planck scale 
Mp = M^Tc- Thus, we have a situation which is inverted with respect to 
normal compactification: four-dimensional gravity at short distances and 
five-dimensional gravity at long distances. 

11 
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2.3 Randall- Sundrum plus brane— induced grav- 
ity 

In this subsection we will investigate what happens when both mechanisms 
are at work simultaneously. The relevant effective action now is 

J dy d^x^iM^R - A) + /" d^^dy6{y) ^/^giM^'r.R - V,) (2.24) 

The crucial observation here is that since the RS solution is ffat on the brane 
, it is not affected by the presence of the localized Einstein term. Thus with 
the RS fine-tuning V^^ = -12KM^ the solution ^T^, ^T^ is still valid. 

Now the equation for the static (scalar) graviton propagator is modified 
to 

M\-e-^^Vl - d'y - ^A'dy - r,6{y)Vl)G{x; y) = 5{y)5'^'\x) (2.25) 

This can be Fourier transformed along the x* coordinated obtaining 

M3(-e-2V -dl~ AAdy + r,5{y)p')G{p- y) = 5{y) (2.26) 

Imposing the symmetry G(p, —y) = G{p,y) we obtain the solution 

G{p,y)=Bw'K,i^) (2.27) 

where the constant B is given by 

We investigate the force mediated by the graviton fluctuations on the 
brane by evaluating 

We distinguish two separated cases. 

For Ktc ^ 1 the potential exhibits four-dimensional behavior ~ 1/r for 
all distances, with an effective Planck scale M^i ~ M^Tc- 

For Ktc ^ 1 we find five-dimensional behavior ~ 1/r^ for the potential 
for distances Tc ^^ r <^ 1/-^) while for r ^ 1/K or r <^ re the potential 
displays four-dimensional behavior. For short distances the Planck scale is 
Mp ~ M^r^i while for long distances the effective Planck scale is M^ 
M^/K. 



r-^ 
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Chapter 3 



Preliminaries 



3.1 The Standard Model of particle physics 
and supersymmetry 

3.1.1 The SM of particle physics 

The Standard Model Lagrangian 

The Standard Model is a gauge theory of the microscopic interactions. 
The strong interaction part is a gauge theory, based on the gauge group 
SU{3) and is described by the Lagrangian 

Lsu, = -\f;,F^'^'' + J2 Qra^ rp q^ (3.1) 

r 

where 

-^/li/ = d^_,Gl - duG''^ - gsfijk G^ G^ (3.2) 

is the field strength tensor for the gluon fields G^, i = 1, ■ ■ ■ , 8, with gs 
the QCD gauge coupling constant, and the structure constants fijk {i,j, k = 
1, ■ ■ • ,8) are defined by 

[\\\']=2if,,kX'' (3.3) 

where the SU{3) A matrices are defined in Table 13.11 The F'^ term leads 
to three and four-point gluon self-interactions. The second term in Lsu3 
is the gauge covariant derivative for the quarks: qr is the r^^ quark fiavor, 
a, P = 1,2,3 are color indices, and 

^M/3 = (D,U = d,6.P + Ws G; L^^ (3.4) 
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Table 3.1: The SU{3) matrices. 

where the quarks transform according to the triplet representation matrices 
L' = X^/2. The color interactions are diagonal in the flavor indices, but in 
general change the quark colors. They are purely vector (parity conserving). 
There are no bare mass terms for the quarks in ()3.H) . These would be allowed 
by QCD alone, but are forbidden by the chiral symmetry of the electroweak 
part of the theory. The quark masses will be generated later by spontaneous 
symmetry breaking. There are in addition effective ghost and gauge-fixing 
terms which enter into the quantization of both the SU{3) and electroweak 
lagrangians, and there is the possibility of adding an (unwanted) term which 
violates CP invariance. 

The electroweak theory is based on the SU{2) x f/(l) gauge group. Its 
Lagrangian is 

LsU2XUi = -^gauge + L^ + Lf + Lyukawa (3-5) 



The gauge part is 



^gauge = -\f;^F'^- - \b,,B^- (3.6) 



where VF^, i = 1, 2, 3 and B^ are respectively the SU{2) and U{1) gauge 
fields, with field strength tensors 

Bfj,u = d^B^ — duB^ 

F,. = d,Wl-dM;-ge,,uWf,W^ (3.7) 

where g{g') is the SU{2) (f/(l)) gauge coupling and eijk is the totally antisym- 
metric symbol. The SU{2) fields have three and four-point self- interactions. 
5 is a f/(l) field associated with the weak hypercharge Y = Q — T3, where 
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Q and T3 are respectively the electric charge operator and the third compo- 
nent of weak SU{2). It has no self-interactions. The B and W^ fields will 
eventually mix to form the photon and Z boson. 

The scalar part of the lagrangian is 

L^ = {D^v)^D^^ - V{^) (3.8) 

where V^ = I ) is a complex Higgs scalar, which is a doublet under 
SU{2) with f/(l) charge Y^p = +|. The gauge covariant derivative is 

D,^ = (^d, + tg^W; + '-^B,^ ^ (3.9) 

where the r* are the Pauli matrices. The square of the covariant derivative 
leads to three and four-point interactions between the gauge and scalar fields. 

V{ip) is the Higgs potential. The combination of SU{2) x U{1) invariance 
and renormalizability restricts V to the form 

V{^) = +/iVV + A(<^V)' (3.10) 

For /i^ < there will be spontaneous symmetry breaking. The A term 
describes a quartic self-interaction between the scalar fields. Vacuum stability 
requires A > 0. 

The fermion term is 

F 

+dl^i]/)dlR + el^i]/)el^) (3.11) 



m=l 



In fl3.11|) m is the family index, F > 3 is the number of families, and L{R) 
refer to the left (right) chiral projections ipiiR) = (1 T 75)V'/2- The left- 
handed quarks and leptons 



7/0 \ / 7/0 



<ilL=[j^)^ Cl = [ l-o ) (3.12) 



L 



transform as SU{2) doublets, while the right-handed fields m^^, (i^^, and 
^mR ^^^ singlets. Their U{1) charges are F^^ = 1,1^^ = — |, F^^ = q^. The 
superscript refers to the weak eigenstates, i.e. fields transforming according 
to definite SU{2) representations. They may be mixtures of mass eigenstates 
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(flavors). The quark color indices a = r, g, b have been suppressed. The 
gauge covariant derivatives are 



D.e^mR = id,-w'B,)el^ (3.13) 

from which one can read off the gauge interactions between the W and B 
and the fermion fields. The different transformations of the L and R fields 
(i.e. the symmetry is chiral) is the origin of parity violation in the elec- 
troweak sector. The chiral symmetry also forbids any bare mass terms for 
the fermions. 

The last term in ()3.5j) is 



-L'Yukawa — 2_^ [^mnlmLf'^mR + ^mnQmLfd^^R + ^mn''mnf^nR\ + H.C. 
m,n=l 

(3.14) 
where the matrices Tmn describe the Yukawa couplings between the single 
Higgs doublet, ip, and the various fiavors m and n of quarks and leptons. 
One needs representations of Higgs fields with Y = ^ and — | to give masses 
to the down quarks, the electrons, and the up quarks. The representation 
(f^^ has y = — |, but transforms as the 2* rather than the 2. However, in 
SU{2) the 2* representation is related to the 2 by a similarity transformation, 

and Lp = ir'^ip^ = j _ I transforms as a 2 with Y"^ = — i. All of the 

masses can therefore be generated with a single Higgs doublet if one makes 
use of both (p and (p. The fact that the fundamental and its conjugate 
are equivalent does not generalize to higher unitary groups. Furthermore, 
in supersymmetric extensions of the standard model supersymmetry forbids 
the use of a single Higgs doublet in both ways in the lagrangian, and one 
must add a second Higgs doublet. Similar statements apply to most theories 
with an additional U{1) gauge factor, i.e. a heavy Z' boson. 
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Spontaneous Symmetry Breaking 

Gauge invariance (and therefore renormalizability) does not allow mass terms 
in the lagrangian for the gauge bosons or for chiral fermions. Massless gauge 
bosons are not acceptable for the weak interactions, which are known to be 
short-ranged. Hence, the gauge invariance must be broken spontaneously, 
which preserves the renormalizability. The idea is simply that the lowest 
energy (vacuum) state does not respect the gauge symmetry and induces 
effective masses for particles propagating through it. 

Let us introduce the complex vector 

v= {0\ip\0) = constant (3.15) 

which has components that are the vacuum expectation values of the various 
complex scalar fields, v is determined by rewriting the Higgs potential as a 
function of v, V{(f) — ^ V{v), and choosing v such that V is minimized. That 
is, we interpret v as the lowest energy solution of the classical equation of 
motion. The quantum theory is obtained by considering fluctuations around 
this classical minimum, ip = v + ^p'. 

The single complex Higgs doublet in the standard model can be rewritten 
in a Hermitian basis as 



where (pi = (p] represent four hermitian flelds. In this new basis the Higgs 
potential becomes 

which is clearly O4 invariant. Without loss of generality we can choose the 
axis in this four-dimensional space so that (0|(/9j|0) = 0, i = 1,2,4 and 
(0|(^3|0) = u. Thus, 

\/(^) ^ \/(t;) = i/xV + Iaz.^ (3.18) 

which must be minimized with respect to u. For /i^ > the minimum occurs 
at z/ = 0. That is, the vacuum is empty space and SU{2) x U{1) is unbroken 
at the minimum. On the other hand, for /i^ < the z/ = symmetric point is 
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V. 
V 



unstable, and the minimum occurs at some nonzero value of v which breaks 
the SU{2) X U{1) symmetry. The point is foun d by requiring 

\/'(zy) = Zy(^2 ^ Az/2) = (3.19) 

which has the solution u = (— /i^/A) at the minimum. (The solution for 
—u can also be transformed into this standard form by an appropriate O4 
transformation.) The dividing point /i^ = cannot be treated classically. It 
is necessary to consider the one loop corrections to the potential, in which 
case it is found that the symmetry is again spontaneously broken. 

We are interested in the case fi^ < 0, for which the Higgs doublet is 

replaced, in first approximation, by its classical value V' "^ 75 I 

The generators L^, L^, and L^ — Y are spontaneously broken (e.g. L^v 7^ 0). 
On the other hand, the vacuum carries no electric charge {Qv = {L^ + Y)v = 
0), so the UiQ of electromagnetism is not broken. Thus, the electroweak 
SU{2) X f/(l) group is spontaneously broken down, SU{2) x Uiy -^ Uiq. 

To quantize around the classical vacuum, write ip = v + ip', where (/?' are 
quantum fields with zero vacuum expectation value. To display the physical 
particle content it is useful to rewrite the four hermitian components of ip' 
in terms of a new set of variables using the Kibble transformation: 

i^ is a hermitian field which will turn out to be the physical Higgs scalar. If 
we had been dealing with a spontaneously broken global symmetry the three 
hermitian fields ^* would be the massless pseudoscalar Goldstone bosons that 
are necessarily associated with broken symmetry generators. However, in a 
gauge theory they disappear from the physical spectrum. To see this it is 
useful to go to the unitary gauge 

,^,-.,-.E..V^^(^;^) (3,21) 

in which the Goldstone bosons disappear. In this gauge, the scalar covariant 
kinetic energy term takes the simple form 



(D.^YD'^^ = i(Oz/) 



|rw; + |5, 



H terms 

V 



Ml 
M^W+^W^ + -^Z^Z^ + H terms (3.22) 
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where the kinetic energy and gauge interaction terms of the physical H par- 
ticle have been omitted. Thus, spontaneous symmetry breaking generates 
mass terms for the W and Z gauge bosons 

v2 
Z = -sm9wB + cos OwW^ (3.23) 

The photon field 

A = cos 9wB + sin 9wW^ (3.24) 

remains massless. The masses are 

Mw = Y (3-25) 

and 

M, 



Mz = V9' + 9"^7^ = ^^ (3.26) 

2 cos Ow 

where the weak angle is defined by tan^vi^ = 9'/ 9- One can think of the 
generation of masses as due to the fact that the W and Z interact constantly 
with the condensate of scalar fields and therefore acquire masses, in analogy 
with a photon propagating through a plasma. The Goldstone boson has 
disappeared from the theory but has reemerged as the longitudinal degree of 
freedom of a massive vector particle. 

It can be shown that Gp/V^ ~ 9^/8M^, where Gp = 1.16639(2) x 
10^^ GeV~^ is the Fermi constant determined by the muon lifetime. The 
weak scale z/ is therefore 

z/ = 2Mw/9 ^ {V2Gf)-^^^ ~ 246 GeV (3.27) 

Similarly, g = e/ sinOw, where e is the electric charge of the positron. Hence, 
to lowest order 

M^ = M, cos e^ r. i-^/V2G,Y/' ^3 28) 

smb'vK 

where a ~ 1/137.036 is the fine structure constant. Using sin^ 9w ~ 0.23 
from neutral current scattering, one expects Mw ~ 78 GeV, and Mz ~ 
89 GeV. (These predictions are increased by ~ (2 — 3) GeV by loop correc- 
tions.) The W and Z were discovered at CERN by two groups (UAl and 
UA2) in 1983. Subsequent measurements of their masses and other proper- 
ties have been in perfect agreement with the standard model expectations 
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(including the higher-order corrections), as is described in the articles of by 
Schaile and Einsweiler. 

After symmetry breaking the Higgs potential becomes 

V{v) = -^^-^^'H' + X^H' + \H^ (3.29) 

The third and fourth terms represent the cubic and quartic interactions of 
the Higgs scalar. The second term represents a (tree-level) mass 



Mh = V-2/i2 = V2Az/ (3.30) 

The weak scale is given in ()3.27|1 . but the quartic Higgs coupling A is un- 
known, so Mh is not predicted. A priori, A could be anywhere in the range 
< A < cxD. There is now an experimental lower limit Mh > 60 GeV from 
LEP. Otherwise, the decay Z ^ Z*H would have been observed (There are 
also theoretical lower limits on Mh in the ~ 10 GeV range, depending on 
rrit, when higher-order corrections are included). 

3.1.2 Basics of Super symmetry 

The motivation for supersymmetry 

It is widely accepted that the Standard Model of gauge interactions describ- 
ing the laws of physics at the weak scale is extraordinarily successful. The 
agreement between theory and experimental data is very good. Yet, we 
believe that the present structure is incomplete. Only to remind a few draw- 
backs, the theory has too many parameters, it does not describe the fermion 
masses and why the number of generations is three. It contains fundamental 
scalars, something difficult to reconcile with our current understanding of 
non- super symmetric field theory. Finally, it does not incorporate gravity. 

It is tempting to speculate that a new (but yet undiscovered) symme- 
try, supersymmetry j^, may provide answers to these fundamental questions. 
Supersymmetry is the almost universally accepted framework for construct- 
ing extensions of the Standard Model. Supersymmetry can be formulated 
either as a global or a local symmetry. In the latter case it includes gravity, 
and is therefore called supergravity. Supersymmetry is the only framework 
in which we seem to be able to understand light fundamental scalars. It 
addresses the question of parameters: first, unification of gauge couplings 
works much better with than without supersymmetry; second, it is easier to 
attack questions such as fermion masses in supersymmetric theories, in part 
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simply due to the presence of fundamental scalars. Supersymmetry seems to 
be intimately connected with gravity. So there are a number of theoretical 
arguments that suggest that nature might be supersymmetric, and that su- 
persymmetry might manifest itself at energies of order the weak interaction 
scale. 

The supersymmetry algebra and supermultiplets 

We begin with some basics, that apply to both global supersymmetry and 
supergravity. 

In the low-energy regime, phenomenology requires the type of supersym- 
metry known as A^ = 1 (one generator). In this section, we present some 
features of A^ = 1 supersymmetric theories, that are likely to be relevant for 
inflation. 

The basic supersymmetry algebra is given by 

{g,,g^} = 2a^p^ (3.31) 

where Qa and Qh are the supersymmetry generators (bars stand for conju- 
gate), a and (3 run from 1 to 2 and denote the two-component Weyl spinors 
(quantities with dotted indices transform under the (0, |) representation of 
the Lorentz group, while those with undotted indices transform under the 
(|,0) conjugate representation), a^^ is a matrix four vector, a^ = (— l,a) 
and P^ is the generator of spacetime displacements (four-momentum). 

The chiral and vector superfields are two irreducible representations of 
the supersymmetry algebra containing fields of spin less than or equal to one. 
Chiral fields contain a Weyl spinor and a complex scalar; vector fields contain 
a Weyl spinor and a (massless) vector. In superspace a chiral superfield may 
be expanded in terms of the Grassmann variable 9 (the fermionic coordinates) 

0(a:, 9) = (j){x) + V29iIj{x) + 9^F{x) (3.32) 

where 9"^ = eai39°'9^. Here x denotes a point in spacetime, 0(a;) is the complex 
scalar, ip the fermion, and F is an auxiliary field. As in this expression, we 
shall generally use the same symbol to represent a superfield and its scalar 
component. Under a supersymmetry transformation with anticommuting 
parameter (, the component fields transform as 

50 = y/2(^ (3.33) 

Sip = V2CF + V2ia^Cdf,(j) (3.34) 

6F = -V2idf,ilja^C (3.35) 

21 



Preliminaries 



Here and in the following, for any generic two-component Weyl spinor A, A 
indicates the complex conjugate of A. For a gauge theory one has to introduce 
vector superfields and the physical content is most transparent in the Wess- 
Zumino gauge. In this gauge and for the simplest case of an abelian group 
f/(l), the vector superfield may be written as 

V = -Oa^'OA^ + ioWx - ie^ex + \e'^e'^D (3.36) 

Here A^ is the gauge field, A^ is the gaugino, and D is an auxiliary field. 
The analog of the gauge invariant field strength is a chiral field: 

W^ = -iK + e^D - f (a'^a"^)„F^, + ^V^S^A'^ (3.37) 

where F^^ = d^A^ — d^A^, and where a^ = (—1, — (?). Regarding the super- 
symmetry transformations, let us just note that 

5\ = KD + C<y^a''F^^ (3.38) 

Global supersymmetry is defined as invariance under these transforma- 
tions with (^ independent of spacetime position, and local supersymmetry 
(supergravity) as invariance with C, depending on spacetime position. In 
the latter case one has to introduce another supermultiplet containing the 
graviton and the gravitino. 

Global supersymmetry may be regarded as a limit of supergravity, in 
which roughly speaking gravity is made negligible by taking Mpi to infinity. 
For most purposes it is a good approximation if the vevs of all relevant scalar 
fields and auxiliary fields are much less than M-p\. 

The Lagrangian of global supersymmetry 

We focus first on global supersymmetry, with the usual restriction that it be 
renormalizable. 

To write down the action for a set of chiral superfields, (pi, transforming 
in some representation of a gauge group G, one introduces, for each gauge 
generator, a vector superfield, V"". Defining the matrix V = T"'Va, where 
T" are the hermitian generators of the gauge group G in the representation 
defined by the scalar fields and excluding the possible Fayet-Iliopoulos term 
to be discussed later, the most general renormalizable lagrangian, written in 
superspace, is then 

£ = ^ /" (fed^9(j)le^(Pn + i / d^^W^ + [ d^dW{(j)n) + h.c. (3.39) 
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where in the adjoint representation Tr(T"T'') = M"* and W{(f)n{x,0)) is a 
fundamental object known as superpotential. The corresponding function of 
the scalar components (pnix), denoted by the same name and symbol, is a 
holomorphic function of the 0„. For simplicity, we shall pretend that there is 
a single gauge f/(l) interaction, with coupling constant g. This is adequate 
since such an interaction is the only one that we consider in detail. (To be 
precise, we consider a U{1) with a Fayet-Iliopoulos term.) In the case of 
several f/(l)'s, there are no cross-terms in the potential from the D-terms, 
i.e. Vd is simply expressed as J2n(^D)n- 

To write this down in terms of component fields, we need the covariant 
derivative 

D^ = d^- '-gA^ (3.40) 

In terms of the component fields, the lagrangian takes the form: 



n 






+ J2^^ 



+ C.C. (3.41) 

At the end of the second line, g„ are the f/(l)-charges of the fields 0„. The 
equations of motion for the auxiliary fields F„ and D are the constraints: 

D = -fE^«l'^"l' (3.43) 

n 

Eq. ()3.4H) contains the gauge invariant kinetic terms for the various 
fields, which specify their gauge interactions. It also contains, after having 
made use of Eqs. ()3.42|1 and ()3.43j) . the scalar field potential, 

V = Vf + Vd (3.44) 

Vf = J^li^nl' (3.45) 

n 

Vd = \d^ (3.46) 
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This separation of the potential into an F term and a D term is crucial 
for inflation model-building, especially when it is generalized to the case of 
supergravity. 

The potential specifies the masses of the scalar fields, and their interac- 
tions with each other. The first term in the third line specifies the interactions 
of gaugino and scalar fields, while the second specifies the masses of the chi- 
ral fermions and their interactions with the scalars. All of these non-gauge 
interactions are called Yukawa couplings. 

To have a renormalizable theory, W is at most cubic in the fields, corre- 
sponding to a potential which is at most quartic. 

From the above expressions, in particular Eq. ()3.45p . one sees that the 
overall phase of W is not physically significant. An internal symmetry can 
either leave W invariant, or alter its phase. The latter case corresponds to 
what is called an R-symmetry. Because W is holomorphic, the internal sym- 
metries restrict its form much more than is the case for the actual potential 
V . In particular, terms in W of the form ^rrKpl or m(f)i(j)2i which would 
generate a mass term m^|0ip in the potential, are usually forbidden.^ As 
a result, scalar particles usually acquire masses only from the vevs of scalar 
fields (i.e. from the spontaneous breaking of an internal symmetry) and from 
supersymmetry breaking. The same applies to the spin-half partners of scalar 
fields, with the former contribution the same in both cases. 

In the case of a U{1) gauge symmetry, one can add to the above la- 
grangian what is called a Fayet-Iliopoulos term 

-2^ f (fed'^e V (3.47) 

This corresponds to adding a contribution —$, to the D field, so that ()3.43j) 
becomes 

^ = -|$^gn|0nr-e (3.48) 

n 

The D term of the potential therefore becomes 

From now on, we shall use a more common notation, where ^ and the 



-•^An exception is the /i term of the MSSM, jiHuHe,, which gives mass to the Higgs 
fields. 
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charges are redefined so that 

^D = ^/f5^gn|0„r + n (3.50) 

This is equivalent to 

D = -giyQn\<Pn\' + a (3.51) 



\ n 



A Fayet-Ihopoulos term may be present in the underlying theory from the 
very beginning, or appears in the effective theory after some heavy degrees 
of freedom have been integrated out. 

Spontaneously broken global susy 

Global supersymmetry breaking may be either spontaneous or explicit. How- 
ever here we shall discuss only the first case. For spontaneous breaking, the 
lagrangian is supersymmetric as given in the last subsection. But the gen- 
erators Qa fail to annihilate the vacuum. Instead, they produce a spin-half 
field, which may be either a chiral field tpa or a gauge field Aq,. The con- 
dition for spontaneous susy breaking is therefore to have a nonzero vacuum 
expectation value for {Qa,'ipf3} or {Qa,^i3}- 

The former quantity is defined by Eq. ()3.34|) , and the latter by Eq. fl3.38|) . 
The quantities d^^cj) and F^,^ contain derivatives of fields, and are supposed 
to vanish in the vacuum. It follows that susy is spontaneously broken if, and 
only if, at least one of the auxiliary fields F„ or D has a non-vanishing vev. 

In the true vacuum, one defines the scale Ms of global supersymmetry 
breaking by 

Ml = J2m' + lD' (3.52) 

n 

or equivalently 

MI = V (3.53) 

(In the simplest case D vanishes and there is just one F„.) 

When we go to supergravity, part of V is still generated by the super- 
symmetry breaking terms, but there is also a contribution — 3|iyp/M|[. This 
allows V to vanish in the true vacuum as is (practically) demanded by ob- 
servation. 
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During inflation, V is positive so the negative term is smaller than the 
susy-breaking terms. In most models of inflation it is negligible. In any case, 
V is at least as big as the susy breaking term, so the search for a model of 
inflation is also a search for a susy-breaking mechanism in the early Universe. 

Spontaneous symmetry breaking can be either tree-level (already present 
in the lagrangian) or dynamical (generated only by quantum effects like con- 
densation). The spontaneous breaking in general breaks the equality between 
the scalar and spin-i masses, in each chiral supermultiplet. But at tree level 
the breaking satisfies a simple relation, which can easily be derived from the 
lagrangian ()H.4H1 . Ignoring mass mixing for simplicity, one finds in the case 
of symmetry breaking by an F-term, 



I]«i 



+ ml,-2ml)=0 (3.54) 



Here n labels the chiral supermultiplets, m„f is the fermion mass while rrini 
and mn2 are the scalar masses. In the case of symmetry breaking by a D 
term, coming from a ^(1), the right hand side of Eq. ()3.54j) becomes DTrQ. 
But in order to cancel gauge anomalies, it is strongly desirable that TrQ = 
which recovers Eq. ()3.54|) . 



3.2 The Standard Model of Cosmology, the 
early Universe and inflation 

3.2.1 The SM of cosmology 

The Robertson- Walker Metric 

Cosmology jl] as the application of general relativity (GR) to the entire 
universe would seem a hopeless endeavor were it not for a remarkable fact - 
the universe is spatially homogeneous and isotropic on the largest scales. 

"Isotropy" is the claim that the universe looks the same in all directions. 
Direct evidence comes from the smoothness of the temperature of the cosmic 
microwave background. "Homogeneity" is the claim that the universe looks 
the same at every point. It is harder to test directly, although some evidence 
comes from number counts of galaxies. More traditionally, we may invoke the 
"Copernican principle," that we do not live in a special place in the universe. 
Then it follows that, since the universe appears isotropic around us, it should 
be isotropic around every point; and a basic theorem of geometry states that 
isotropy around every point implies homogeneity. 
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We may therefore approximate the universe as a spatially homogeneous 
and isotropic three-dimensional space which may expand (or, in principle, 
contract) as a function of time. The metric on such a spacetime is necessarily 
of the Robert son- Walker (RW) form. 

Therefore, the most general spacetime metric consistent with homogene- 
ity and isotropy is 



ds^ = -de + a^{t) [dp" + f{p) [de^ + sin^ edcj)'')] 
where the three possibilities for /(p) are 

/(p) = {sin(p), p, sinh(p)} 



(3.55) 



(3.56) 



This is a purely geometric fact, independent of the details of general relativ- 
ity. We have used spherical polar coordinates (p, 6*, 0), since spatial isotropy 
implies spherical symmetry about every point. The time coordinate t, which 
is the proper time as measured by a comoving observer (one at constant spa- 
tial coordinates), is referred to as cosmic time, and the function a{t) is called 
the scale factor. 

There are two other useful forms for the RW metric. First, a simple 
change of variables in the radial coordinate yields 



ds^ 



where 



-dt^ + a^{t) 



+1 



-1 



dr"^ 



kr^ 



+ r^ {de^ + sin^ ^rf0^) 



if /(p) = sin(p) 

if /(P)=P 

if /(p) = sinh(p) 



(3.57) 



(3.58) 



Geometrically, k describes the curvature of the three-dimensional space. 
k = +1 corresponds to positively curved spatial sections (locally isometric 
to 3-spheres); k = corresponds to local flatness, and k = —1 corresponds 
to negatively curved (locally hyperbolic) spatial sections. 

Note that we have not chosen a normalization such that ao = 1. We 
are not free to do this and to simultaneously normalize \k\ = 1, without 
including explicit factors of the current scale factor in the metric. In the flat 
case, where /c = 0, we can safely choose oq = 1. 

A second change of variables, which may be applied to either ()3.55|) 
or ()3.57|) . is to transform to conformal time, r, via 



it) 



dt' 

a(¥) 



(3.59) 
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Applying this to ()3.57|) yields 
ds^ = aHr 



dr^ + ^^ + r^ (d9' + sin^ 9d^') 
1 — kr^ ^ ' 



(3.60) 



where we have written aij^) = a[t{T)] as is conventional. The confornial time 
does not measure the proper time for any particular observer, but it does 
simplify some calculations. 

A particularly useful quantity to define from the scale factor is the Hubble 
parameter (sometimes called the Hubble constant), given by 

H = - (3.61) 

a 

The Hubble parameter relates how fast the most distant galaxies are receding 
from us to their distance from us via Hubble's law, 

V ^ Hd. (3.62) 

This is the relationship that was discovered observationally by Edwin Hubble, 
and has been verified to high accuracy by modern observational methods. 



The Friedmann Equations 

As mentioned, the RW metric is a purely kinematic consequence of requiring 
homogeneity and isotropy of our spatial sections. We next turn to dynamics, 
in the form of differential equations governing the evolution of the scale factor 
a(t). These will come from applying Einstein's equation, 

Ri,u - -^Rgfiu = 8itGT^^ (3.63) 

to the RW metric. 

Before diving right in, it is useful to consider the types of energy-mome- 
ntum tensors T^i, we will typically encounter in cosmology. For simplicity, 
and because it is consistent with much we have observed about the universe, 
it is often useful to adopt the perfect fluid form for the energy-momentum 
tensor of cosmological matter. This form is 

Tf,u = {p + p)U^U,+pg^, (3.64) 

where If^ is the fluid four-velocity, p is the energy density in the rest frame of 
the fluid and p is the pressure in that same frame. The pressure is necessarily 
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isotropic, for consistency with the RW metric. Similarly, fluid elements will 
be comoving in the cosmological rest frame, so that the normalized four- 
velocity in the coordinates of ()3.57|) will be 

[/'^ = (1, 0, 0, 0) (3.65) 

The energy-momentum tensor thus takes the form 

where gij represents the spatial metric (including the factor of a^). 

Armed with this simplified description for matter, we are now ready to 
apply Einstein's equation ()3.63|) to cosmology. Using ()3.57|1 and ()3.64|1 . one 
obtains two equations. The first is known as the Friedmann equation, 

where an overdot denotes a derivative with respect to cosmic time t and i 
indexes all different possible types of energy in the universe. This equation 
is a constraint equation, in the sense that we are not allowed to freely specify 
the time derivative d; it is determined in terms of the energy density and 
curvature. The second equation, which is an evolution equation, is 



a + Ka) =-^'^Ep-5;^ (3-68) 



It is often useful to combine ()3.67p and ()3.68|) to obtain the acceleration 
equation 

- = ^5Z(p. + 3p,) (3.69) 

i 

In fact, if we know the magnitudes and evolutions of the different energy 
density components pi, the Friedmann equation ()3.67|) is sufficient to solve 
for the evolution uniquely. The acceleration equation is conceptually useful, 
but rarely invoked in calculations. 

The Friedmann equation relates the rate of increase of the scale factor, as 
encoded by the Hubble parameter, to the total energy density of all matter 
in the universe. We may use the Friedmann equation to define, at any given 
time, a critical energy density, 

Pc^^ (3.70) 
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for which the spatial sections must be precisely flat (fc = 0). We then define 
the density parameter 

^total = — , (3-71) 

Pc 

which allows us to relate the total energy density in the universe to its local 
geometry via 

^total > 1 ^ k = +l 

aotai = l ^ /e = (3.72) 

^total < 1 ^ k= -I 

It is often convenient to define the fractions of the critical energy density in 
each different component by 

a = - (3.73) 

Pc 

Energy conservation is expressed in GR by the vanishing of the covariant 
divergence of the energy-momentum tensor, 

V^T"" = (3.74) 

Applying this to our assumptions - the RW metric ()3.57j) and perfect-fluid 
energy-momentum tensor ()3.(i4|l - yields a single energy-conservation equa- 
tion, 

p + 3H{p + p)=Q (3.75) 

This equation is actually not independent of the Friedmann and acceleration 
equations, but is required for consistency. It implies that the expansion of the 
universe (as specified by H) can lead to local changes in the energy density. 
Note that there is no notion of conservation of "total energy," as energy can 
be interchanged between matter and the spacetime geometry. 

One final piece of information is required before we can think about 
solving our cosmological equations: how the pressure and energy density are 
related to each other. Within the fluid approximation used here, we may 
assume that the pressure is a single-valued function of the energy density 
P = P{p)- It is often convenient to define an equation of state parameter, w, 
by 

p = wp (3.76) 

This should be thought of as the instantaneous definition of the parameter w; 
it need represent the full equation of state, which would be required to cal- 
culate the behavior of fluctuations. Nevertheless, many useful cosmological 
matter sources do obey this relation with a constant value of w. For example, 
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w = corresponds to pressureless matter, or dust - any collection of massive 
non-relativistic particles would qualify. Similarly, w = 1/3 corresponds to 
a gas of radiation, whether it be actual photons or other highly relativistic 
species. 

A constant w leads to a great simplification in solving our equations. In 
particular, using ()3.75j) . we see that the energy density evolves with the scale 
factor according to 

/^(«) ^ ^(^ (3.77) 

Note that the behaviors of dust {w = 0) and radiation {w = 1/3) are consis- 
tent with what we would have obtained by more heuristic reasoning. Consider 
a fixed comoving volume of the universe - i.e. a volume specified by fixed 
values of the coordinates, from which one may obtain the physical volume at 
a given time t by multiplying by a{t)^. Given a fixed number of dust particles 
(of mass m) within this comoving volume, the energy density will then scale 
just as the physical volume, i.e. as a(t)~^, in agreement with ()3.77|) . with 
w = 0. 

To make a similar argument for radiation, first note that the expansion 
of the universe (the increase of a{t) with time) results in a shift to longer 
wavelength A, or a redshift, of photons propagating in this background. A 
photon emitted with wavelength Ag at a time te, at which the scale factor 
is Qe = a{te) is observed today {t = to, with scale factor oq = a(to)) at 
wavelength Aq, obeying 

^ = ^ = 1 + ;^ (3.78) 

Ae 0,e 

The redshift z is often used in place of the scale factor. Because of the 
redshift, the energy density in a fixed number of photons in a fixed comoving 
volume drops with the physical volume (as for dust) and by an extra factor 
of the scale factor as the expansion of the universe stretches the wavelengths 
of light. Thus, the energy density of radiation will scale as a(t)~^, once again 
in agreement with ()3.77|1 . with w = 1/3. 

Thus far, we have not included a cosmological constant A in the gravi- 
tational equations. This is because it is equivalent to treat any cosmological 
constant as a component of the energy density in the universe. In fact, adding 
a cosmological constant A to Einstein's equation is equivalent to including 
an energy-momentum tensor of the form 

T,. = -^9,. (3.79) 
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This is simply a perfect fluid with energy-momentum tensor ()3.64j) with 

A 



PA = -PA (3.80) 

so that the equation-of-state parameter is 

wa = -1 (3.81) 

This implies that the energy density is constant, 

Pa = constant (3.82) 

Thus, this energy is constant throughout spacetime; we say that the cosmo- 
logical constant is equivalent to vacuum energy. 

Similarly, it is sometimes useful to think of any nonzero spatial curvature 
as yet another component of the cosmological energy budget, obeying 



P: 



3k 



curv 



SirGa^ 
Pcu. = ^ (3.83) 

so that 

Wcnr. = -1/3 (3.84) 

It is not an energy density, of course; Pcurv is simply a convenient way to keep 
track of how much energy density is lacking, in comparison to a flat universe. 



Flat Universes 

It is much easier to find exact solutions to cosmological equations of motion 
when k = 0. Fortunately for us, nowadays we are able to appeal to more than 
mathematical simplicity to make this choice. Indeed, modern cosmological 
observations, in particular precision measurements of the cosmic microwave 
background, show the universe today to be extremely spatially flat. 

In the case of flat spatial sections and a constant equation of state pa- 
rameter w, we may exactly solve the Friedmann equation ()3.77p to obtain 

/ , N 2/3(l+t«) 

a{t)=aoi-] (3.85) 
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Type of Energy 


p{a) 


ait) 


Dust 

Radiation 

Cosmological Constant 


constant 





Table 3.2: A summary of the behaviors of the most important sources of 
energy density in cosmology. The behavior of the scale factor applies to 
the case of a fiat universe; the behavior of the energy densities is perfectly 
general. 

where Qq is the scale factor today, unless w = —1, in which case one obtains 
a{t) oc e^*. Applying this result to some of our favorite energy density sources 
yields Table EISl 

Note that the matter- and radiation-dominated fiat universes begin with 
a = 0; this is a singularity, known as the Big Bang. We can easily calculate 
the age of such a universe: 



I 







da 



,0 aH{a) 3{l + w)Ho 
Unless w is close to —1, it is often useful to approximate this answer by 



(3.86) 



t, 







H, 



-1 





(3.87) 



It is for this reason that the quantity Hq^ is k;nown as the Hubble time, and 
provides a useful estimate of the time scale for which the universe has been 
around. 

Horizons 

One of the most crucial concepts to master about FRW models is the ex- 
istence of horizons. This concept will prove useful in understanding the 
shortcomings of what we are terming the standard cosmology. 

Suppose an emitter, e, sends a light signal to an observer, o, who is at 
r = 0. Setting 6 = constant and (p = constant and working in conformal 
time, for such radial null rays we have To — t = r. In particular this means 
that 

To-Te = re (3.88) 

Now suppose Tg is bounded below by f^; for example, fg might represent the 
Big Bang singularity. Then there exists a maximum distance to which the 
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observer can see, known as the particle horizon distance, given by 

rph{ro) = T-o - fe (3.89) 

Similarly, suppose Tq is bounded above by Tq. Then there exists a limit 
to spacetime events which can be influenced by the emitter. This limit is 
known as the event horizon distance, given by 

reh{ro) = fo-Te (3.90) 

These horizon distances may be converted to proper horizon distances at 
cosmic time t, for example 

/■* dt' 

dn = a{T)rpi^ = a{T){T -Te) =a{t) — — (3.91) 

Just as the Hubble time Hq^ provides a rough guide for the age of the 
universe, the Hubble distance cHq^ provides a rough estimate of the horizon 
distance in a matter- or radiation-dominated universe. 

3.2.2 The early Universe 

In this subsection we use what we know of the laws of physics and the universe 
today to infer conditions in the early universe. Early times were character- 
ized by very high temperatures and densities, with many particle species kept 
in (approximate) thermal equilibrium by rapid interactions. We will there- 
fore have to move beyond a simple description of non-interacting "matter" 
and "radiation," and discuss how thermodynamics works in an expanding 
universe. 

Describing Matter 

We have discussed how to describe matter as a perfect fluid, described by an 
energy-momentum tensor 

T^,u = {p + p)U^U,+pg^, (3.92) 

where If^ is the fluid four-velocity, p is the energy density in the rest frame 
of the fluid and p is the pressure in that same frame. The energy-momentum 
tensor is covariantly conserved, 

V^T^^ = (3.93) 

34 



3.2. The Standard Model of Cosmology, the early Universe and inflation 

In a more complete description, a fluid will be characterized by quantities 
in addition to the energy density and pressure. Many fluids have a conserved 
quantity associated with them and so we will also introduce a number flux 
density N^^, which is also conserved 

V^A^'^ = (3.94) 

For non-tachyonic matter N^^ is a timelike 4-vector and therefore we may 
decompose it as 

A^^ = nW (3.95) 

We can also introduce an entropy flux density S'^. This quantity is not 
conserved, but rather obeys a covariant version of the second law of thermo- 
dynamics 

^^,S>' > (3.96) 

Not all phenomena are successfully described in terms of such a local entropy 
vector (e.g. black holes); fortunately, it suffices for a wide variety of fluids 
relevant to cosmology. 

The conservation law for the energy-momentum tensor yields, most im- 
portantly, equation ()3.75p . which can be thought of as the first law of ther- 
modynamics 

dU = TdS - pdV (3.97) 

with dS = 0. 

It is useful to resolve S^ into components parallel and perpendicular to 
the fiuid 4-velocity 

S" = sU" + s^ (3.98) 

where Sfj^lf^ = 0. The scalar s is the rest-frame entropy density which, up to 
an additive constant (that we can consistently set to zero), can be written 
as 

s = ^ (3.99) 

In addition to all these quantities, we must specify an equation of state, 
and we typically do this in such a way as to treat n and s as independent 
variables. 

Particles in Equilibrium 

The various particles inhabiting the early universe can be usefully charac- 
terized according to three criteria: in equilibrium vs. out of equilibrium 
(decoupled), bosonic vs. fermionic, and relativistic (velocities near c) vs. 
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non-relativistic. In this section we consider species which are in equihbrium 
with the surrounding thermal bath. 

Let us begin by discussing the conditions under which a particle species 
will be in equilibrium with the surrounding thermal plasma. A given species 
remains in thermal equilibrium as long as its interaction rate is larger than 
the expansion rate of the universe. Roughly speaking, equilibrium requires it 
to be possible for the products of a given reaction to have the opportunity to 
recombine in the reverse reaction and if the expansion of the universe is rapid 
enough this won't happen. A particle species for which the interaction rates 
have fallen below the expansion rate of the universe is said to have frozen 
out or decoupled. If the interaction rate of some particle with the background 
plasma is F, it will be decoupled whenever 

T<^H (3.100) 

where the Hubble constant H sets the cosmological timescale. 

As a good rule of thumb, the expansion rate in the early universe is 
"slow," and particles tend to be in thermal equilibrium (unless they are very 
weakly coupled). This can be seen from the Friedmann equation when the 
energy density is dominated by a plasma with p ~ T^; we then have 

Thus, the Hubble parameter is suppressed with respect to the temperature 
by a factor of T/Mpi. At extremely early times (near the Planck era, for 
example), the universe may be expanding so quickly that no species are 
in equilibrium; as the expansion rate slows, equilibrium becomes possible. 
However, the interaction rate F for a particle with cross-section a is typically 
of the form 

T = n{av) , (3.102) 

where n is the number density and v a typical particle velocity. Since n oc a~^, 
the density of particles will eventually dip so low that equilibrium can once 
again no longer be maintained. In our current universe, no species are in 
equilibrium with the background plasma (represented by the CMB photons). 

Now let us focus on particles in equilibrium. For a gas of weakly- 
interacting particles, we can describe the state in terms of a distribution 
function /(p), where the three- momentum p satisfies 

E^{p)=m^ + \p\^ (3.103) 
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Relativistic 


Relativistic 


Non-relativistic 




Bosons 


Fermions 


(Either) 


Hi 


'if9.T^ 


Uj ^2 3i^ 


9^{^f"e-'/^ 


Pi 


i>9^T^ 


(1\ ila-T^ 


miTii 


Pi 


\pi 


Ipi 


riiT < Pi 



Table 3.3: Number density, energy density, and pressure, for species in ther- 
mal equilibrium. 

The distribution function characterizes the density of particles in a given 
momentum bin. (In general it will also be a function of the spatial position 
X, but we suppress that here.) The number density, energy density, and 
pressure of some species labeled i are given by 



Hi 
pi 
Pi 



9i 



(27r)3 

9i 
(27r)3 

9i 



Mp)d^P 
Eip)f,ip)d'p 



iPl 



(27r)3y 3E{p) 



Mp)d'p 



(3.104) 



where Qi is the number of spin states of the particles. For massless photons we 
have g^ = 2, while for a massive vector boson such as the Z we have 9z = '^■ 
In the usual accounting, particles and antiparticles are treated as separate 
species; thus, for spin- 1/2 electrons and positrons we have ge- = ge+ = 2. 
In thermal equilibrium at a temperature T the particles will be in either 
Fermi-Dirac or Bose-Einstein distributions, 



/(P) 



1 



eE{p)/T ± 1 



(3.105) 



where the plus sign is for fermions and the minus sign for bosons. 

We can do the integrals over the distribution functions in two oppo- 
site limits: particles which are highly relativistic (T ^ m) or highly non- 
relativistic {T <^ m). The results are shown in Table ESI in which ( is the 
Riemann zeta function, and ({3) ~ 1.202. 

From Table 13.31 we can extract several pieces of relevant information. 
Relativistic particles, whether bosons or fermions, remain in approximately 
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equal abundances in equilibrium. Once they become non-relativistic, how- 
ever, their abundance plummets, and becomes exponentially suppressed with 
respect to the relativistic species. This is simply because it becomes pro- 
gressively harder for massive particle-antiparticle pairs to be produced in a 
plasma with T <^m. 

It is interesting to note that, although matter is much more dominant 
than radiation in the universe today, since their energy densities scale differ- 
ently the early universe was radiation-dominated. We can write the ratio of 
the density parameters in matter and radiation as 



^ = ^f-')=^(l + .r (3.106) 

"R "Ro \ao/ "RO 

The redshift of matter-radiation equality is thus 

l + ^cq = ^~3xl03 (3.107) 

"RO 

This expression assumes that the particles that are non-relativistic today 
were also non-relativistic at Zeq; this should be a safe assumption, with the 
possible exception of massive neutrinos, which make a minority contribution 
to the total density. 

At this point we should stress that even decoupled photons maintain a 
thermal distribution; this is not because they are in equilibrium, but simply 
because the distribution function redshifts into a similar distribution with a 
lower temperature proportional to 1/a. We can therefore speak of the "effec- 
tive temperature" of a relativistic species that freezes out at a temperature 
Tf and scale factor a/: 



Tr'ia) = T, 



e) p- 



For example, neutrinos decouple at a temperature around 1 ~ MeV; shortly 
thereafter, electrons and positrons annihilate into photons, dumping energy 
(and entropy) into the plasma but leaving the neutrinos unaffected. Con- 
sequently, we expect a neutrino background in the current universe with a 
temperature of approximately 2K, while the photon temperature is 3K. 

A similar effect occurs for particles which are non-relativistic at decou- 
pling, with one important difference. For non-relativistic particles the tem- 
perature is proportional to the kinetic energy |mf^, which redshifts as 1/a^. 
We therefore have 
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In either case we are imagining that the species freezes out while relativistic 
/non-relativistic and stays that way afterward; if it freezes out while rela- 
tivistic and subsequently becomes non-relativistic, the distribution function 
will be distorted away from a thermal spectrum. 

The notion of an effective temperature allows us to define a corresponding 
notion of an effective number of relativistic degrees of freedom, which in 
turn permits a compact expression for the total relativistic energy density. 
The effective number of relativistic degrees of freedom (as far as energy is 
concerned) can be defined as 

9.= Y.4f)+l E 3.(f) (3.110) 

bosons ^ fermions ^ ^ 

(The temperature T is the actual temperature of the background plasma, 
assumed to be in equilibrium.) Then the total energy density in all relativistic 
species comes from adding the contributions of each species, to obtain the 
simple formula 

We can do the same thing for the entropy density. From ()3.99|) . the entropy 
density in relativistic particles goes as T^ rather than T'^, so we define the 
effective number of relativistic degrees of freedom for entropy as 

9*s=Y.3^{Jf)U E 3.{Jf) (3.112) 

bosons fermions ^ 

The entropy density in relativistic species is then 

s = %9*sT'' (3.113) 

Numerically, g^, and g^:S will typically be very close to each other. In the 
Standard Model, we have 

C 100, T > 300 MeV 
9* ~ ^*s ~ < 10, 1 MeV < T < 300 MeV (3.114) 

[3, T < 1 MeV 

The events that change the effective number of relativistic degrees of free- 
dom are the QCD phase transition at 300 ~ MeV, and the annihilation of 
electron/positron pairs at 1 ~ MeV. 

Because of the release of energy into the background plasma when species 
annihilate, it is only an approximation to say that the temperature goes as 
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T oc 1/a. A better approximation is to say that the comoving entropy density 
is conserved, 

socfl-^ (3.115) 

This will hold under all forms of adiabatic evolution; entropy will only be pro- 
duced at a process like a first-order phase transition or an out-of-equilibrium 
decay. (In fact, we expect that the entropy production from such processes 
is very small compared to the total entropy, and adiabatic evolution is an 
excellent approximation for almost the entire early universe.) Combining 
entropy conservation with the expression ()3.113|) for the entropy density in 
relativistic species, we obtain a better expression for the evolution of the 
temperature. 

The temperature will consistently decrease under adiabatic evolution in an 
expanding universe, but it decreases more slowly when the effective number 
of relativistic degrees of freedom is diminished. 



Toc(7;/Vi (3.116) 



Thermal Relics 

As we have mentioned, particles typically do not stay in equilibrium forever; 
eventually the density becomes so low that interactions become infrequent, 
and the particles freeze out. Since essentially all of the particles in our current 
universe fall into this category, it is important to study the relic abundance 
of decoupled species. (Of course it is also possible to obtain a significant relic 
abundance for particles which were never in thermal equilibrium; examples 
might include baryons produced by GUT baryogenesis, or axions produced 
by vacuum misalignment.) In this subsection we will typically neglect factors 
of order unity. 

We have seen that relativistic, or hot, particles have a number density 
that is proportional to T^ in equilibrium. Thus, a species X that freezes out 
while still relativistic will have a number density at freeze-out Ty given by 

nxiTf)r~^Tf (3.117) 

Since this is comparable to the number density of photons at that time, 
and after freeze-out both photons and our species X just have their number 
densities dilute by a factor a(t)~^ as the universe expands, it is simple to 
see that the abundance of X particles today should be comparable to the 
abundance of CMB photons, 

nxo ~ «70 ~ 10^ cm"^ (3.118) 
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We express this number as 10^ rather than 411 since the roughness of our 
estimate does not warrant such misleading precision. The leading correction 
to this value is typically due to the production of additional photons subse- 
quent to the decoupling of X; in the Standard Model, the number density of 
photons increases by a factor of approximately 100 between the electro weak 
phase transition and today, and a species which decouples during this period 
will be diluted by a factor of between 1 and 100 depending on precisely when 
it freezes out. So, for example, neutrinos which are light {m,y < MeV) have 
a number density today of n,, = 115 cm~^ per species, and a corresponding 
contribution to the density parameter (if they are nevertheless heavy enough 
to be nonrelativistic today) of 

no. ^ (5^) .- (3.119) 

(In this final expression we have secretly taken account of the missing numer- 
ical factors, so this is a reliable answer.) Thus, a neutrino with m„ ~ 10~^ eV 
would contribute ^2,^ ~ 2 x 10~^. This is large enough to be interesting with- 
out being large enough to make neutrinos be the dark matter. That's good 
news, since the large velocities of neutrinos make them free-stream out of 
overdense regions, diminishing primordial perturbations and leaving us with 
a universe which has much less structure on small scales than we actually 
observe. 

Now consider instead a species X which is nonrelativistic or cold at the 
time of decoupling. It is much harder to accurately calculate the relic abun- 
dance of a cold relic than a hot one, simply because the equilibrium abun- 
dance of a nonrelativistic species is changing rapidly with respect to the 
background plasma, and we have to be quite precise following the freeze- 
out process to obtain a reliable answer. The accurate calculation typically 
involves numerical integration of the Boltzmann equation for a network of 
interacting particle species; here, we cut to the chase and simply provide a 
reasonable approximate expression. If o"o is the annihilation cross-section of 
the species X at a temperature T = rrix, the final number density in terms 
of the photon density works out to be 

nx{T < Tf) ^-^^n^ (3.120) 

Since the particles are nonrelativistic when they decouple, they will certainly 
be nonrelativistic today, and their energy density is 

Px =mxnx (3.121) 
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We can plug in numbers for the Hubble parameter and photon density to 
obtain the density parameter, 

Numerically, when h = c = 1 we have 1 GeV~ 2 x 10^^^ cm, so the photon 
density today is n^ ~ 100 cm^^ ~ 10^'^^ GeV^^. The Hubble constant is 
Hq ~ 10~^^ GeV, and the Planck mass is Mpi ~ lO"*^^ GeV, so we obtain 

^x --^ 97 (3.123) 

o-o(109 GeV^) ^ ' 

It is interesting to note that this final expression is independent of the mass 
mx of our relic, and only depends on the annihilation cross-section; that's 
because more massive particles will have a lower relic abundance. Of course, 
this depends on how we choose to characterize our theory; we may use vari- 
ables in which ctq is a function of mx, in which case it is reasonable to say 
that the density parameter does depend on the mass. 

One candidate for Cold Dark Matter (CDM) is a Weakly Interacting 
Massive Particle (WIMP). The annihilation cross-section of these particles, 
since they are weakly interacting, should be o"o ~ c^w^F' where aw is the 
weak coupling constant and Gp is the the Fermi constant. Using Gp ~ 
(300GeV)-2 and aw ~ 10-^ we get 

do ~ a^G^ ~ 10-9 GeV~^ (3.124) 

Thus, the density parameter in such particles would be 

fix ~ 1 (3.125) 

In other words, a stable particle with a weak interaction cross section nat- 
urally produces a relic density of order the critical density today, and so 
provides a perfect candidate for cold dark matter. A paradigmatic example 
is provided by the lightest supersymmetric partner (LSP), if it is stable and 
supersymmetry is broken at the weak scale. Such a possibility is of great 
interest to both particle physicists and cosmologists, since it may be possi- 
ble to produce and detect such particles in colliders and to directly detect a 
WIMP background in cryogenic detectors in underground laboratories; this 
will be a major experimental effort over the next few years. 

Baryogenesis 

The symmetry between particles and antiparticles, firmly established in col- 
lider physics, naturally leads to the question of why the observed universe is 
composed almost entirely of matter and no primordial antimatter. 
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If large domains of matter and antimatter exist, then annihilations would 
take place at the interface between them. If the typical size of such a domain 
was small enough, then the energy released by these annihilations would 
result in a diffuse 7-ray background and a distortion of the cosmic microwave 
radiation, neither of which is observed. 

While the above considerations put an experimental upper bound on the 
amount of antimatter in the universe, strict quantitative estimates of the 
relative abundances of baryonic matter and antimatter may also be obtained 
from the standard cosmology. The baryon number density does not remain 
constant during the evolution of the universe, instead scaling like a~^, where 
a is the cosmological scale factor. It is therefore convenient to define the 
baryon asymmetry of the universe in terms of the quantity 

r] = — (3.126) 

s 

where s is the entropy density and ub is the difference between the baryon 
number density and the anti-baryon number density. The range of r] consis- 
tent with the observational data is 

2.6 X 10"^° < ?7 < 6.2 X 10"^° (3.127) 

Thus the natural question arises: As the universe cooled from early times to 
today, what processes, both particle physics and cosmological, were respon- 
sible for the generation of this very specific baryon asymmetry? 

As pointed out by Sakharov, a small baryon asymmetry r] may have been 
produced in the early universe if three necessary conditions are satisfied 

• Baryon number (B) violation 

• Violation of C (charge conjugation symmetry) and CP (the combina- 
tion of C and parity) 

• departure from thermal equilibrium 

The first condition should be clear since, starting from a baryon symmet- 
ric universe with 77 = 0, baryon number violation must take place in order 
to evolve into a universe in which 77 does not vanish. The second Sakharov 
criterion is required because, if C and CP are exact symmetries, one can 
prove that the total rate for any process which produces an excess of baryons 
is equal to the rate of the complementary process which produces an ex- 
cess of antibaryons and so no net baryon number can be created. That is 
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to say that the thermal average of the baryon number operator B, which 
is odd under both C and CP, is zero unless those discrete symmetries are 
violated. CP violation is present either if there are complex phases in the La- 
grangian which cannot be reabsorbed by field redefinitions (explicit breaking) 
or if some Higgs scalar field acquires a VEV which is not real (spontaneous 
breaking). We will discuss this in detail shortly. 

Finally, to explain the third criterion, one can calculate the equilibrium 
average of S at a temperature T = 1/(3: 

{B)t = TT{e~^^B) = Ti[{CPT){CPTyh~^^B)] 

= Tr {e-'^"{CPT)-^B{CPT)] = -Tr (e-f^^B) (3.128) 

where we have used that the Hamiltonian H commutes with CPT. Thus 
{B)t = in equilibrium and there is no generation of net baryon number. 

Of the three Sakharov conditions, baryon number violation and C and 
CP violation may be investigated only within a given particle physics model, 
while the third condition - the departure from thermal equilibrium - may be 
discussed in a more general way, as we shall see. Let us discuss the Sakharov 
criteria in more detail. 

Baryon Number Violation 

Grand Unified Theories (GUTs) 5J describe the fundamental interactions by 
means of a unique gauge group G which contains the Standard Model (SM) 
gauge group SU{?>)c ® SU{2)l ® U{1)y. The fundamental idea of GUTs 
is that at energies higher than a certain energy threshold Mqut the group 
symmetry is G and that, at lower energies, the symmetry is broken down to 
the SM gauge symmetry, possibly through a chain of symmetry breakings. 
The main motivation for this scenario is that, at least in supersymmetric 
models, the (running) gauge couplings of the SM unify at the scale Mqut — 
2 X 10^^ GeV, hinting at the presence of a GUT involving a higher symmetry 
with a single gauge coupling. 

Baryon number violation seems very natural in GUTs. Indeed, a general 
property of these theories is that the same representation of G may contain 
both quarks and leptons, and therefore it is possible for scalar and gauge 
bosons to mediate gauge interactions among fermions having different baryon 
number. 

It is well-known that the most general renormalizable Lagrangian invari- 
ant under the SM gauge group and containing only color singlet Higgs fields 
is automatically invariant under global abelian symmetries which may be 
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identified with the baryonic and leptonic symmetries. These, therefore, are 
accidental symmetries and as a result it is not possible to violate B and L 
at tree-level or at any order of perturbation theory. Nevertheless, in many 
cases the perturbative expansion does not describe all the dynamics of the 
theory and, indeed, in 1976 't Hooft realized that nonperturbative effects 
(instantons) may give rise to processes which violate the combination B + L, 
but not the orthogonal combination B — L. The probability of these pro- 
cesses occurring today is exponentially suppressed and probably irrelevant. 
However, in more extreme situations - like the primordial universe at very 
high temperatures - baryon and lepton number violating processes may be 
fast enough to play a significant role in baryogenesis. 

CP violation 

CP violation in GUTs arises in loop-diagram corrections to baryon number 
violating bosonic decays. Since it is necessary that the particles in the loop 
also undergo 5-violating decays, the relevant particles are the X, Y, and H^ 
bosons in the case of SU{5). 

In the electroweak theory things are somewhat different. Since only the 
left-handed fermions are SU{2)l gauge coupled, C is maximally broken in 
the SM. Moreover, CP is known not to be an exact symmetry of the weak 
interactions. This is seen experimentally in the neutral kaon system through 
Ko, Kq mixing. Thus, CP violation is a natural feature of the standard 
electroweak model. 

While this is encouraging for baryogenesis, it turns out that this partic- 
ular source of CP violation is not strong enough. The relevant effects are 
parameterized by a dimensionless constant which is no larger than 10~^°. 
This appears to be much too small to account for the observed BAU and, 
thus far, attempts to utilize this source of CP violation for electroweak baryo- 
genesis have been unsuccessful. In light of this, it is usual to extend the SM 
in some fashion that increases the amount of CP violation in the theory while 
not leading to results that conflict with current experimental data. One con- 
crete example of a well-motivated extension in the minimal supersymmetric 
standard model (MSSM). 

In some scenarios, such as GUT baryogenesis, the third Sakharov con- 
dition is satisfied due to the presence of superheavy decaying particles in a 
rapidly expanding universe. These generically fall under the name of out-of- 
equilibrium decay mechanisms. 

The underlying idea is fairly simple. If the decay rate Tx of the su- 
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per heavy particles X at the time they become nonrelativistic (i.e. at the 
temperature T ~ Mx) is much smaller than the expansion rate of the uni- 
verse, then the X particles cannot decay on the time scale of the expansion 
and so they remain as abundant as photons for T < Mx- In other words, 
at some temperature T > Mx, the superheavy particles X are so weakly 
interacting that they cannot catch up with the expansion of the universe and 
they decouple from the thermal bath while they are still relativistic, so that 
nx ^ riy r^ T^ at the time of decoupling. 

Therefore, at temperature T ~ Mx, they populate the universe with an 
abundance which is much larger than the equilibrium one. This overabun- 
dance is precisely the departure from thermal equilibrium needed to produce 
a final nonvanishing baryon asymmetry when the heavy states X undergo B 
and CP violating decays. 

The out-of-equilibrium condition requires very heavy states: Mx > (10^^ 
-10^^) GeV and Mx > (10^° - 10^^) GeV, for gauge and scalar bosons, 
respectively, if these heavy particles decay through renormalizable operators. 

Since the linear combination B — L is left unchanged by sphaleron tran- 
sitions, the baryon asymmetry may be generated from a lepton asymmetry. 
Indeed, sphaleron transition will reprocess any lepton asymmetry and con- 
vert (a fraction of) it into baryon number. This is because B + L must be 
vanishing and the final baryon asymmetry results to be i? ~ —L. 

In the SM as well as in its unified extension based on the group SU{5), 
B — L is conserved and no asymmetry in B — L can be generated. However, 
adding right-handed Majorana neutrinos to the SM breaks B — L and the 
primordial lepton asymmetry may be generated by the out-of-equilibrium 
decay of heavy right-handed Majorana neutrinos iV£ (in the supersymmetric 
version, heavy scalar neutrino decays are also relevant for leptogenesis). This 
simple extension of the SM can be embedded into GUTs with gauge groups 
containing 5*0(10). Heavy right-handed Majorana neutrinos can also explain 
the smallness of the light neutrino masses via the see-saw mechanism. 



3.2.3 Inflation 

So far we have described what is known as the standard cosmology. This 
framework is a towering achievement, describing to great accuracy the phys- 
ical processes leading to the present day universe. However, there remain 
outstanding issues in cosmology. Many of these come under the heading of 
initial condition problems and require a more complete description of the 
sources of energy density in the universe. The most severe of these problems 
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Fixed Point 


1 + 3w; > 


1 + 3u7 < 


Q = 

n = 1 

fi = oo 


attractor 

repeller 

attractor 


repeller 

attractor 

repeller 



Table 3.4: Behavior of the density parameter near fixed points. 

eventually led to a radical new picture of the physics of the early universe - 
cosmological inflation, which is the subject of this subsection. 

We will begin by describing some of the problems of the standard cos- 
mology. 



The Flatness Problem 

The Friedmann equation may be written as 

k 



fi-1 



H^a 



2n2 



(3.129) 



where for brevity we are now writing Q instead of fitotai- Differentiating this 
with respect to the scale factor, this implies 



da a 



(3.130) 



This equation is easily solved, but its most general properties are all that 
we shall need and they are qualitatively different depending on the sign of 
1 + 3w. There are three fixed points of this differential equation, as given in 
Table m 

Observationally we know that 1] ^ 1 today - i.e. we are very close to 
the repeller of this differential equation for a universe dominated by ordinary 
matter and radiation {w > —1/3). Even if we only took account of the lu- 
minous matter in the universe, we would clearly live in a universe that was 
far from the attractor points of the equation. It is already quite puzzling 
that the universe has not reached one of its attractor points, given that the 
universe has evolved for such a long time. However, we may be more quan- 
titative about this. If the only matter in the universe is radiation and dust, 
then in order to have Q in the range observed today requires (conservatively) 



< 1 - fi < 10 



-60 



(3.131) 
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This remarkable degree of fine tuning is the flatness problem. Within the 
context of the standard cosmology there is no known explanation of this 
fine-tuning. 

The Horizon Problem 

The horizon problem stems from the existence of particle horizons in FRW 
cosmologies, as discussed in a previous subsection. Horizons exist because 
there is only a finite amount of time since the Big Bang singularity, and thus 
only a finite distance that photons can travel within the age of the universe. 
Consider a photon moving along a radial trajectory in a flat universe (the 
generalization to non-flat universes is straightforward). In a flat universe, we 
can normalize the scale factor to 

ao = 1 (3.132) 

without loss of generality. A radial null path obeys 

= ds^ = -dt^ + a^dr'^ (3.133) 

so the comoving (coordinate) distance traveled by such a photon between 
times ti and t2 is 

Ar= r^T (3-134) 

To get the physical distance as it would be measured by an observer at any 
time t, simply multiply by a{t). For simplicity let's imagine we are in a 
matter-dominated universe, for which 



(e 



\ 2/3 

(3.135) 



The Hubble parameter is therefore given by 

2 

-t 

3 



H = 2r' 



= a-^/^i/o (3.136) 

Then the photon travels a comoving distance 

Ar = 2H^^ (v/^ - vV) (3.137) 

The comoving horizon size when a = a* is the distance a photon travels since 
the Big Bang, 

rhor(a*) = 2H^^^, (3.138) 
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The physical horizon size, as measured on the spatial hypersurface at a*, is 
therefore simply 

dhor{a*) = a*rhor(a*) = 2//^^ (3.139) 

Indeed, for any nearly-flat universe containing a mixture of matter and radi- 
ation, at any one epoch we will have 

4or(a*) ~ H;' (3.140) 

where H~^ is the Hubble distance at that particular epoch. This approximate 
equality leads to a strong temptation to use the terms "horizon distance" and 
"Hubble distance" interchangeably; this temptation should be resisted, since 
inflation can render the former much larger than the latter, as we will soon 
demonstrate. 

The horizon problem is simply the fact that the CMB is isotropic to 
a high degree of precision, even though widely separated points on the last 
scattering surface are completely outside each others' horizons. When we look 
at the CMB we were observing the universe at a scale factor ocmb ~ 1/1200; 
meanwhile, the comoving distance between a point on the CMB and an 
observer on Earth is 

Ar = 2H^' {I - ,/E^^) 

f« 2Hq^ (3.141) 

However, the comoving horizon distance for such a point is 

^hor(aCMB) = 2ijQ y/acMB 

^ 6 X IQ-^Hq^ (3.142) 

Hence, if we observe two widely-separated parts of the CMB, they will have 
non-overlapping horizons; distinct patches of the CMB sky were causally 
disconnected at recombination. Nevertheless, they are observed to be at the 
same temperature to high precision. The question then is, how did they know 
ahead of time to coordinate their evolution in the right way, even though they 
were never in causal contact? We must somehow modify the causal structure 
of the conventional FRW cosmology. 

Unwanted Relics 

We have already talked about grand unifled theories (GUTs). If grand uni- 
flcation occurs with a simple gauge group G, any spontaneous breaking of G 
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satisfies tt2{G/H) = 7ri(if) for any simple subgroup H. In particular, break- 
ing down to the standard model will lead to magnetic monopoles 0, since 

7T2{G/H) = TTi{[SU{3) X SU{2) X U{1)]/Zg) = Z (3.143) 

(The gauge group of the standard model is, strictly speaking, [SU{3) x 
SU{2) X U{1)]/Zq. The Zq factor only affects the global structure of the 
group and not the Lie algebra, and thus is usually ignored by particle physi- 
cists.) 

Using the Kibble mechanism, the expected relic abundance of monopoles 
works out to be 

"»~-i°"(i^)'(i^^) P"*' 

This is far too big; the monopole abundance in GUTs is a serious problem 
for cosmology if GUTs have anything to do with reality. 

In addition to monopoles, there may be other model-dependent relics 
predicted by our favorite theory. If these are incompatible with current 
limits, it is necessary to find some way to dilute their density in the early 
universe. 

The General Idea of Inflation 

The horizon problem especially is an extremely serious problem for the stan- 
dard cosmology because at its heart is simply causality. Any solution to this 
problem is therefore almost certain to require an important modification to 
how information can propagate in the early universe. Cosmological inflation 
is such a mechanism. 

Before getting into the details of inflation we will just sketch the general 
idea here. The fundamental idea is that the universe undergoes a period of 
accelerated expansion, defined as a period when d > 0, at early times. The 
effect of this acceleration is to quickly expand a small region of space to a 
huge size, diminishing spatial curvature in the process, making the universe 
extremely close to flat. In addition, the horizon size is greatly increased, 
so that distant points on the CMB actually are in causal contact and un- 
wanted relics are tremendously diluted, solving the monopole problem. As 
an unexpected bonus, quantum fluctuations make it impossible for inflation 
to smooth out the universe with perfect precision, so there is a spectrum of 
remnant density perturbations; this spectrum turns out to be approximately 
scale-free, in good agreement with observations of our current universe. 
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Slowly-Rolling Scalar Fields 

If inflation is to solve the problems of the standard cosmology, then it must be 
active at extremely early times. Thus, we would like to address the earliest 
times in the universe amenable to a classical description. We expect this 
to be at or around the Planck time tp and since Planckian quantities arise 
often in inflation we will retain values of the Planck mass in the equations 
of this section. There are many models of inflation, but because of time 
constraints we will concentrate almost exclusively on the chaotic inflation 
model of Linde. We have borrowed heavily in places here from the excellent 
text of Liddle and Lyth. 

Consider modeling matter in the early universe by a real scalar fleld 0, 
with potential V{(f)). The energy-momentum tensor for is 



Tf,u = (y^,(t>)(yu(t>) - Qt^u 



^^"^(V«0)(V^</') + ^(</>) 



(3.145) 



For simplicity we will specialize to the homogeneous case, in which all quan- 
tities depend only on cosmological time t and set A; = 0. A homogeneous real 
scalar field behaves as a perfect fluid with 

P0 = \^'' + V{<P) (3.146) 

P4> = 1<P'-V{<P) (3.147) 

The equation of motion for the scalar fleld is given by 

'^ + 3-0+^ = (3.148) 

a d(j) 

which can be thought of as the usual equation of motion for a scalar fleld 
in Minkowski space, but with a friction term due to the expansion of the 
universe. The Friedmann equation with such a fleld as the sole energy source 
is 



H 



2 



SttG 



U' + V^(0) 



(3.149) 



A very speciflc way in which accelerated expansion can occur is if the uni- 
verse is dominated by an energy component that approximates a cosmological 
constant. In that case the associated expansion rate will be exponential, as 
we have already seen. Scalar flelds can accomplish this in an interesting way. 
From ()3.14()|1 it is clear that if cfp' <^ V{(j)) then the potential energy of the 
scalar fleld is the dominant contribution to both the energy density and the 
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pressure, and the resulting equation of state is p ~ — p, approximately that of 
a cosnaological constant, the resulting expansion is certainly accelerating. In 
a loose sense, this negligible kinetic energy is equivalent to the fields slowly 
rolling down its potential; an approximation which we will now make more 
formal. 

Technically, the slow-roll approximation for inflation involves neglecting 
the term in ()3.148|) and neglecting the kinetic energy of compared to 
the potential energy. The scalar field equation of motion and the Friedmann 
equation then become 

*^-^ (3.150) 

H'- c. ^V(4,) (3.151) 

where a prime denotes a derivative with respect to 0. 

These conditions will hold if the two slow-roll conditions are satisfied. 
These are 

|e| < 1 

\v\ < 1 (3.152) 



where the slow-roll parameters are given by 



Ml fV'V 
e=^ 77 (3-153) 



y , 

and 

ri^Ml,^ (3.154) 

It is easy to see that the slow roll conditions yield inflation. Recall that 
inflation is defined by a/ a > 0. We can write 



a 



H + H' (3.155) 



so that inflation occurs if 



But in slow-roll 



-^ > -1 (3.156) 

§-,^-e (3.157) 

which will be small. Smallness of the other parameter rj helps to ensure that 
inflation will continue for a sufficient period. 
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It is useful to have a general expression to describe how much inflation 
occurs, once it has begun. This is typically quantified by the number of 
e-folds, defined by 

mt) . ... (^) (3,158, 

Usually we are interested in how many efolds occur between a given field value 
and the field value at the end of inflation 0end, defined by e(0end) = 1- We 
also would like to express A^ in terms of the potential. Fortunately this is 
simple to do via 

^w - "> {^) - r ^ * - 4 c ^ * '"''' 

The issue of initial conditions for inflation is one that is quite subtle and 
we will not get into a discussion of that here. Instead we will remain focused 
on chaotic inflation, in which we assume that the early universe emerges from 
the Planck epoch with the scalar field taking different values in different parts 
of the universe, with typically Planckian energies. There will then be some 
probability for inflation to begin in some places, and we shall focus on those. 



Vacuum Fluctuations and Perturbations 

Recall that the structures - clusters and superclusters of galaxies - we see on 
the largest scales in the universe today, and hence the observed fluctuations 
in the CMB, form from the gravitational instability of initial perturbations 
in the matter density. The origin of these initial fluctuations is an important 
question of modern cosmology. 

Inflation provides us with a fascinating solution to this problem - in 
a nutshell, quantum fluctuations in the inflaton field during the inflationary 
epoch are stretched by inflation and ultimately become classical fluctuations. 
Let's sketch how this works. 

Since inflation dilutes away all matter fields, soon after its onset the uni- 
verse is in a pure vacuum state. If we simplify to the case of exponential 
inflation, this vacuum state is described by the Gibbons-Hawking tempera- 
ture 

H VV , , 

where we have used the Friedmann equation. Because of this temperature, 
the inflaton experiences fluctuations that are the same for each wavelength 
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(50fc = ^GH- Now, these fluctuations can be related to those in the density 
by 

Sp = ^S<p (3.161) 

Inflation therefore produces density perturbations on every scale. The 
amplitude of the perturbations is nearly equal at each wavenumber, but there 
will be slight deviations due to the gradual change in V as the inflaton rolls. 
We can characterize the fluctuations in terms of their spectrum As{k), related 
to the potential via 

Alik) 



V^ 



(3.162) 

k=aH 



where k = aH indicates that the quantity V^/(y'y is to be evaluated at 
the moment when the physical scale of the perturbation A = a/k is equal 
to the Hubble radius H~^. Note that the actual normalization of ()3.162|) is 
convent ion- dependent, and should drop out of any physical answer. 

The spectrum is given the subscript "S" because it describes scalar fluc- 
tuations in the metric. These are tied to the energy-momentum distribution, 
and the density fluctuations produced by inflation are adiabatic — fluctua- 
tions in the density of all species are correlated. The fluctuations are also 
Gaussian, in the sense that the phases of the Fourier modes describing fluc- 
tuations at different scales are uncorrelated. These aspects of inflationary 
perturbations — a nearly scale-free spectrum of adiabatic density fluctua- 
tions with a Gaussian distribution — are all consistent with current observa- 
tions of the CMB and large-scale structure, and have been conflrmed to new 
precision by WMAP and other CMB measurements. 

It is not only the nearly-massless inflaton that is excited during inflation, 
but any nearly-massless particle. The other important example is the gravi- 
ton, which corresponds to tensor perturbations in the metric (propagating 
excitations of the gravitational fleld). Tensor fluctuations have a spectrum 



^^'*' ~ 4 



(3.163) 

k=aH 



The existence of tensor perturbations is a crucial prediction of inflation which 
may in principle be veriflable through observations of the polarization of 
the CMB. Although CMB polarization has already been detected, this is 
only the E-mode polarization induced by density perturbations; the 5-mode 
polarization induced by gravitational waves is expected to be at a much 
lower level, and represents a signiflcant observational challenge for the years 
to come. 
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For purposes of understanding observations, it is useful to parameterize 
the perturbation spectra in terms of observable quantities. We therefore 
write 

Al{k) oc rs-i (3.164) 

and 

Al{k) oc P^ (3.165) 

where ns and nx are the "spectral indices" . They are related to the slow-roll 
parameters of the potential by 

ns = l-6e + 2r] (3.166) 

and 

riT = -2e (3.167) 

Since the spectral indices are in principle observable, we can hope through 
relations such as these to glean some information about the infiaton potential 
itself. 

Our current knowledge of the amplitude of the perturbations already 
gives us important information about the energy scale of inflation. Note 
that the tensor perturbations depend on V alone (not its derivatives), so 
observations of tensor modes yields direct knowledge of the energy scale. If 
large-scale CMB anisotropics have an appreciable tensor component (pos- 
sible, although unlikely), we can instantly derive inflation ~ (10^^ GeV)'^. 
(Here, the value of V being constrained is that which was responsible for 
creating the observed fluctuations; namely, 60 e-folds before the end of infla- 
tion.) This is remarkably reminiscent of the grand uniflcation scale, which 
is very encouraging. Even in the more likely case that the perturbations 
observed in the CMB are scalar in nature, we can still write 

^iiLon - e^/'lO^' GeV (3.168) 

where e is the slow-roll parameter deflned in ()3.153|) . Although we expect 
e to be small, the 1/4 in the exponent means that the dependence on e is 
quite weak; unless this parameter is extraordinarily tiny, it is very likely that 
KiLon ~ 10^^-10^^ GeV. 

We should note that this discussion has been phrased in terms of the 
simplest models of inflation, featuring a single canonical, slowly-rolling scalar 
fleld. A number of more complex models have been suggested, allowing for 
departures from the relations between the slow-roll parameters and observ- 
able quantities; some of these include hybrid inflation, inflation with novel 
kinetic terms, the curvaton model, low-scale models, brane inflation and 
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models where perturbations arise from modulated coupling constants. This 
list is necessarily incomplete, and continued exploration of the varieties of 
inflationary cosmology will be a major theme of theoretical cosmology into 
the foreseeable future. 

Reheating and Preheating 

Clearly, one of the great strengths of inflation is its ability to redshift away 
all unwanted relics, such as topological defects. However, inflation is not 
discerning, and in doing so any trace of radiation or dust-like matter is sim- 
ilarly redshifted away to nothing. Thus, at the end of inflation the universe 
contains nothing but the inflationary scalar field condensate. How then does 
that matter of which we are made arise? How does the hot big bang phase 
of the universe commence? How is the universe reheated? 

Inflation ends when the slow-roll conditions are violated and, in most 
models, the field begins to fall towards the minimum of its potential. Ini- 
tially, all energy density is in the inflaton, but this is now damped by two 
possible terms. First, the expansion of the universe naturally damps the 
energy density. More importantly, the inflaton may decay into other parti- 
cles, such as radiation or massive particles, both fermionic and bosonic. To 
take account of this one introduces a phenomenological decay term F^ into 
the scalar field equation. The inflaton undergoes damped oscillations and 
decays into radiation which equilibrates rapidly at a temperature known as 
the reheat temperature Trh- 
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Cosmological evolution, inflation and dark mat- 
ter in brane cosmology 



57 



Chapter 4 



Axino dark matter 



We discuss dark matter in the brane world scenario. We work in the Randall- 
Sundrum type II brane world and assume that the lightest supersymmetric 
particle is the axino. We find that the axinos can play the role of cold dark 
matter in the universe, provided that the five- dimensional Planck mass is 
bounded both from below and from above. This is possible for higher reheat- 
ing temperatures compared to the conventional four-dimensional cosmology 
due to a novel expansion law for the universe. 



4.1 Introduction 

There are good theoretical reasons for which particle physics proposes that 
new exotic particles must exist. The most compelling solution to the strong 
CP-problem of quantum chromodynamics (QCD), which can be stated as 
"why is the parameter in QCD so small?", is the one proposed by Peccei 
and Quinn [2j. An additional global, chiral symmetry is introduced, now 
known as Peccei-Quinn (PQ) symmetry, which is spontaneously broken at 
the PQ scale /„ > 10® GeV |H]. Since U{1)pq is a spontaneously broken 
global symmetry, there must be a Nambu-Goldstone boson associated with 
this symmetry. However, because U{1)pq suffers from a chiral anomaly, 
this boson is not massless but acquires a small mass. The pseudo Nambu- 
Goldstone boson associated with this spontaneous symmetry breaking is the 
axion 9J, which has not yet been detected. On the other hand, supersymme- 
try (SUSY) is an ingredient that appears in many theories for physics beyond 
the standard model. SUSY solves the hierarchy problem and predicts that 
every particle we know should be escorted by its superpartner. The axino is 
the superpartner of the axion. In order for the supersymmetric solution of the 
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hierarchy problem to work, it is necessary that the SUSY becomes manifest 
at relatively low energies, less than a few TeV, and therefore the required 
superpartners must have masses below this scale (for supersymmetry and 
supergravity see e.g. [ID,)- 

One of the theoretical problems in modern cosmology is to understand 
the nature of cold dark matter in the universe. There are good reasons, 
both observational and theoretical, to suspect that a fraction of 0.22 of the 
energy density in the universe is in some unknown "dark" form ^Tj. Many 
lines of reasoning suggest that the dark matter consists of some new, as 
yet undiscovered, massive particle which experiences neither electromagnetic 
nor color interactions. In SUSY models which are realized with R-parity 
conservation the lightest supersymmetric particle (LSP) is stable. A popular 
cold dark matter candidate is the LSP, provided that it is electrically and 
color neutral. Certainly the most theoretically developed LSP is the lightest 
neutralino. However, there are other dark matter candidates as well, for 
example the gravitino and the axino. In this work we assume that the axino 
is the LSP. Axinos are special because they have unique properties: They 
are very weekly interacting and their mass can span a wide range, from very 
small (~ eV) to large (~ GeV) values. What is worth stressing is that, in 
contrast to those of the neutralino and the gravitino, axino mass does not 
have to be of the order of the SUSY breaking scale in the visible sector, 
MsusY ~ lOOGeV — ITeV. The first paper to show that the axinos can be 
CDM was ^2] • There are however some early works on axino cosmology (see 

e.g. nni, mi). 

We believe that some time in its early history, the universe experienced 
an inflationary phase ^1- According to the inflationary paradigm, during 
the slow-roll phase of inflation the universe undergoes a rapid expansion, and 
consequently any initial population of axinos is diluted away. After slow-roll 
a reheating phase follows and leads the universe to the radiation era of the 
standard hot Big-Bang cosmology of temperature Tr. As the PQ symmetry 
is restored at fa, we consider only values of T^ up to the PQ scale, which we 
take to be fa = W^^ GeV. Another important scale is the temperature To at 
which axinos decouple from the thermal bath. For Tr > To, there has been 
an early phase in which axinos were in thermal equilibrium with the thermal 
bath. The axino density parameter is then given by the equilibrium number 
density [12] 

If we require that Qah'^ ~ 0.1 then the axino mass rria ~ 0.2 keV. For 
an axino mass in the range rria < 1 keV, 1 keV < rria < 100 keV and 
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rria > lOOkeV, we refer to hot, warm and cold axino dark matter respectively. 
So we see that for Tr> Td, axinos can only be hot dark matter. For Tr < T^, 
the axinos are out of thermal equilibrium so that the production mechanisms 
have to be considered in detail. 

In order to generate a large enough abundance of axinos, one needs to 
repopulate the universe (after inflation) with them. There are two generic 
ways of achieving this. First, they can be generated through thermal produc- 
tion (TP), namely via scattering and decay processes of ordinary particles 
and sparticles still in the thermal bath. Second, axinos may also be produced 
via non-thermal production mechanisms (NTP) possibly present during the 
reheating phase. In ^\ the authors considered both NTP and TP and they 
found that Tr had to be relatively low, below some 10^ GeV. However NTP 
mechanisms are strongly model dependent and we shall not consider them 
here. In ^E] the authors using specific techniques (the hard thermal loop re- 
summation technique pH] together with the Braaten-Yuan prescription [T9] ) 
computed the thermal production rate of axinos in supersymmetric QCD 
and evaluated the relic axino abundance. They found that axinos provide 
the density of cold dark matter observed by WMAP for relatively small re- 
heating temperature after inflation Tr < 10^ GeV, essentially in agreement 
with [17j. Such a low reheating temperature excludes some models for in- 
flation and the baryon asymmetry in the universe has to be explained by a 
mechanism that works efficiently at relatively small temperatures, excluding 
thermal leptogenesis ^20j (for thermal leptogenesis in brane world cosmology 
see EH)- 

The purpose of this paper is to show that this fact can be dealt with in the 
context of the brane world scenario. Our brane world model is the Randall- 
Sundrum type II model (RSII) J22j, and in fact its supersymmetric extended 
model |22j. However, the cosmological solution of this extended model is the 
same as that in the non-supersymmetric model, since Einstein's equations 
belong to the bosonic part. The RSII model offers a novel expansion law for 
the observable four-dimensional universe. We find that the axino abundance 
today is proportional to the transition temperature, at which the modified 
expansion law in the brane world cosmology connects to the standard one, 
rather than the reheating temperature after inflation as in the standard cos- 
mology. This means that even though the reheating temperature can be very 
high, the axinos can play the dominant part of the cold dark matter in the 
universe. Other works that discuss dark matter in brane cosmology are j21] . 

Let us see in more detail the thermal production of axinos (in standard 
cosmology). We assume that after inflation axinos are far from thermal 
equilibrium. With the axino number density ria being much smaller than the 
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photon number density n^, the evolution of Uq, with cosmic time t can be 
described by the Boltzmann equation 

^ + 3Hn^ = C^ (4.2) 

at 

where C^ is the coUision term, while the second term on the left-hand side 
accounts for the dilution of the axinos due to the expansion of the universe 
described by the Hubble parameter H . It is convenient to define the dimen- 
sionless quantity 

1^^ = - (4.3) 

s 

where s is the entropy density for the relativistic degrees of freedom in the 

primordial plasma 

27r2 
s{T) = Kjj{T)—T^ (4.4) 

with heff(T) ~ QeffiT) in the radiation dominated epoch and Qejf counts 
the total number of effectively massless degrees of freedom (those species 
with mass mi ^ T). When all the degrees of freedom are relativistic 
Qeff = 915/4 = 228.75. Replacing the cosmic time t with the temperature 
T, the number density ria with the number-to-entropy ratio 1^ and using 
conservation of the entropy per comoving volume (see for example the first 
reference in JS]), the Boltzmann equation can be cast into the form 

(^^a ^ Ca{T) 

dT ~ Ts{T)H{T) ^ ' 

where H{T) is the Hubble parameter as a function of the temperature T for 
the radiation dominated era 



HiT) = JtMlL_ (4.6) 

^ ^ \ 90 Mpi ^ ^ 

where Mpi = 2.4 x 10^^ GeV is the reduced Planck mass. In terms of the 
number-to-entropy ratio Y^, the axino density parameter is given by 

f^, 2 _ P&h^ _ mgngh'^ _ maYas{To)h'^ 

Per Per Per 

Here we make use of the following values 

To = 2.73K = 2.35 x 10"^^ GeV (4.8) 

KffiTo) = 3.91 (4.9) 

p^Jh^ = 8.1 X IQ-^'^ GeV^ (4.10) 
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The collision term Ca has been computed in supersymmetric QCD by the 
authors of [121 

where A'^c = 3, nj = 6 and g is the QCD coupling constant 

1 3 , / T ^ ^ ~^'^ 



The collision term Ca assumes that gluons and gluinos are in thermal equi- 
librium and hence the expression that gives Ca is only valid for T > rrig. The 
thermal axino production proceeds basically during the hot radiation domi- 
nated epoch, that is at temperatures above that at matter-radiation equality 
Teg. Integrating the Boltzmann equation the axino yield at the present tem- 
perature of the universe To is given by 

or 

3C(3)45V90(iV,2 _ 1)^6 (^0.4336^; + (iV, + nj) In (^-^)) M^^i 

"^ °^ 2/..,,V^4096vrio/2 

(4.14) 
and finally the axino density parameter is obtained 

"" = ^'^^ '« [—) (oTgS?) [wg^) (^— j ("^' 

Considering /„ = 10^^ GeV , axinos can be cold dark matter for masses 
rria > 100 keV and reheating temperatures Tr < 10^ CeV [TE] . 

Recently the brane world models have been attracting a lot of attention 
as a novel higher dimensional theory. In these models, it is assumed that the 
standard model particles are confined on a 3-brane while gravity resides in the 
whole higher dimensional spacetime. The model first proposed by Randall 
and Sundrum (RSII) [22j, is a simple and interesting one, and its cosmological 
evolutions have been intensively investigated |23I2H1I2ZI- According to that 
model, our four-dimensional universe is realized on the 3-brane with a positive 
tension located at the UV boundary of five- dimensional AdS spacetime. In 
the bulk there is just a cosmological constant A5, whereas on the brane 
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there is matter with energy-moment um tensor r^,^. Also, the five-dimensional 
Planck mass is denoted by M5 and the brane tension is denoted by T. If 
Einstein's equations hold for the five- dimensional bulk, then it has been 
shown in |2H] that the effective four-dimensional Einstein's equations induced 
on the brane can be written as 

G^u + Kg^iu = -^Tfiu + (xFa)^^^!^ ~ ^f^^ (4.16) 



m; ■ -Ml 



where g^u is the induced metric on the brane, 7r^,y = 12 ''" ''"mi^ + 8 9^^'' '^"P '^ ~ 
^Tfia^l^ — ^t"^ dfif, A4 is the effective four-dimensional cosmological constant, 
rup is the usual four-dimensional Planck mass and E^^, = C^^^ n^ n^ g^ gl is 
a projection of the five-dimensional Weyl tensor Ca/spa, where n" is the unit 
vector normal to the brane. The tensors tt^^ and E^^, describe the infiuence 
of the bulk in brane dynamics. The five- dimensional quantities are related 
to the corresponding four-dimensional ones through the relations 

mp = A^—Mi (4.17) 

and 

In a cosmological model in which the induced metric on the brane gfj_u has 
the form of a spatially flat Friedmann- Robert son- Walker model, with scale 
factor a{t), the Friedmann-like equation on the brane has the generalized 
form (see e.g. the second reference in ^^) 

where C is an integration constant arising from Efj_^. The cosmological con- 
stant term and the term linear in p are familiar from the four-dimensional 
conventional cosmology. The extra terms, i.e the "dark radiation" term and 
the term quadratic in p, are there because of the presence of the extra di- 
mension. Adopting the Randall- Sundrum fine-tuning 

As = ^ (4.20) 

^ 6M| ^ ^ 

the four-dimensional cosmological constant vanishes. So the generalized 
Friedmann equation takes the final form 

3 V PoJ a 
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where 

Po = QQttGMI (4.22) 

with G the Newton's constant. The second term proportional to p^ and 
the dark radiation are new ingredients in the brane world cosmology and 
lead to a non-standard expansion law. The dark radiation term is severely 
constrained by the success of the Big Bang Nucleosynthesis (BBN), since 
the term behaves like an additional radiation at the BBN era [22] • So, for 
simplicity, we neglect the term in the following analysis. The five-dimensional 
Planck mass is also constrained by the BBN, which is roughly estimated as 
M5 > lOTeV [SO]- A more stringent constraint may be obtained by requiring 
that relative corrections to the Newtonian law of gravity should be small on 
scales r > 1 mm. This gives M5 > 10^ GeV jHT] . 

One can see that the evolution of the early universe can be divided into 
two eras. In the low energy regime p <C Po the first term dominates and we 
recover the usual Friedmann equation of the conventional four- dimensional 
cosmology. In the high density regime po ^ p the second term dominates 
and we get an unconventional expansion law for the universe. In between 
there is a transition temperature Tj for which p{Tt) = pQ. The transition 
temperature Tf is determined as 

r. = 1.6 X 10' f MVV ''^V^GeV (4.23) 



geffj yW^GeV 

once M5 is given. Using the transition temperature the generalized Fried- 
mann-like equation (for the radiation era) can be rewritten in the form 

(4.24) 

with Hst the Hubble parameter in standard four-dimensional Big-Bang cos- 
mology. Assuming a transition temperature Tr 3> Tt and Tt ^ Tgg, the 
following integral can be computed to a very good approximation 

Tr 1 pTt 1 pTji 1 

dT ^=^ = dT ^=^ + dT -^=^ ~ 2Tt (4.25) 




Therefore, the axino yield resulting on integrating the Boltzmann equation 
in brane cosmology is 

^-(^°) - £ "^ r.(r;£!(r) ^ - T.^!^,iT., '^' (-«' 
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The expression is only valid for Tt > nig. So we see that essentially the 
reheating temperature Tr in the axino parameter density Q^ is replaced by 
the transition temperature Tt. Therefore, axinos can play the role of cold 
dark matter in the universe for 

m^ > 100 keV, Tj < 2 x 10^ GeV (4.27) 

independently of the reheating temperature. This is the main result in this 
work. At this point we would like to stress the fact that the axino abundance 
does depend on the reheating temperature, but only through the coupling 
constant g{T^). The function g{T) is a very slow-varying function of the 
temperature (for example, g(T = 10^ GeV) = 0.986 and g(T = 10^° GeV) = 
0.852), so practically we can consider reheating temperatures of the order of 
~ 10^° GeV and neglect the dependence on it. Note that according to the 
analysis of |S2I; with a transition temperature Tt < 10^ GeV the gravitino 
problem can be avoided. It is interesting to note that in order that the axinos 
can play the role of cold dark matter in the universe, the five-dimensional 
Planck mass M5 can only take values in a range between an upper limit and 
a lower limit. If M5 becomes too high, the transition temperature Tj will 
be higher than ~ 10^ GeV and this sets an upper bound for M5: M5 < 
2.9 X 10^° GeV. On the other hand, if M5 becomes too low, the reheating 
temperature will be smaller than the gluino mass, nig ~ 1 TeV. Thus, we 
get a lower bound for M5: M5 > 1.8 x 10^ GeV. So we find a window for the 
five-dimensional Planck mass 

1.8 X 10^ GeV < M5 < 2.9 x 10^° GeV (4.28) 

It is interesting to note that this range for the five-dimensional Planck mass 
is compatible with the bounds mentioned before coming from the BBN and 
modifications to Newton's law, namely that M5 > 10*^ GeV. 

To illustrate the above ideas let us present a specific example. We con- 
sider the case in which the five- dimensional Planck mass is M5 = 10^° GeV. 
Then the transition temperature Tt is found to be Tt = 4 x 10^ GeV. We 
also assume that the reheating temperature is Tr = 10^° GeV. If the axinos 
are to be the cold dark matter in the universe, their parameter density has 
to be Qah"^ = 0.113. This happens for rua — 511 keV. Our treatment is valid 
as long as axinos are never in thermal equilibrium after inflation. One can 
easily check that this is the case. For that we need to compare the Hubble 
parameter H{T) to the rate T{T) of the reaction that maintain the axinos in 
thermal equilibrium. For T < f^ the reaction rate is |14j 

^,3 



Tfofe^ («9) 



' a 
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We can see that after inflation H{T) ^ r(T) for all the values of M5 in the 
allowed range. 

Summarizing, we have discussed dark matter in Randall-Sundrum type 
11 brane world assuming that the axino is the LSP. We have seen that axinos 
can be the dominant part of the cold dark matter in the universe if their 
mass rria > 100 keV and the transition temperature Tt < 2 x 10® GeV 
independently of the reheating temperature Tr after inflation (provided that 
Tr ^ Tt, which is true for Tfi ~ 10^° GeV). Therefore, in contrast to 
the case for the conventional four- dimensional cosmology, high values for Tr 
such as 10^*^ GeV are allowed, in accord with most inflationary models and 
baryogenesis through leptogenesis. 
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Sneutrino inflation 



We discuss sneutrino inflation in the brane-world scenario. We work in the 
Randall- Sundrum type II brane-world, generalized with the introduction of 
the Gauss-Bonnet (GB) term, a correction to the effective action in string 
theories. We find that a viable inflationary model is obtained with a reheating 
temperature appropriate to lead to the right baryon asymmetry and render 
the gravitino safe for cosmology. In speciflc realizations we satisfy all the 
observational constaints without the unnaturally small Yukawa couplings 
required in other related approaches. 



5.1 Introduction 

Inflation jT^j has become the standard paradigm for the early Universe, be- 
cause it solves some outstanding problems present in the standard Hot Big- 
Bang cosmology, like the flatness and horizon problems, the problem of un- 
wanted relics, such as magnetic monopoles, and produces the cosmological 
fluctuations for the formation of the structure that we observe today. The re- 
cent spectacular CMB data from the WMAP satellite [TT| I^ have strengthen 
the inflationary idea, since the observations indicate an almost scale-free 
spectrum of Gaussian adiabatic density fluctuations, just as predicted by 
simple models of inflation. According to chaotic inflation with a potential 
for the inflaton fleld of the form V = (l/2)m^0^, the WMAP normaliza- 
tion condition requires for the inflaton mass m that m = 1.8 x 10^^ GeV [51| . 
However, a yet unsolved problem about inflation is that we do not know how 
to integrate it with ideas in particle physics. For example, we would like to 
identify the inflaton, the scalar fleld that drives inflation, with one of the 
known flelds of particle physics. 
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One of the most exciting experimental results in the last years has been 
the discovery of neutrino oscillations j33|- These results are nicely explained if 
neutrinos have a small but finite mass jSEj- The simplest models of neutrino 
masses invoke heavy gauge-singlet neutrinos that give masses to the light 
neutrinos via the seesaw mechanism [2Z]. If we require that light neutrino 
masses ~ 10~^ to 10"'^ eV, as indicated by the neutrino oscillations data, we 
find that the heavy singlet neutrinos weight ~ 10^° to 10^^ GeV [SHI, a range 
that includes the value of the infiaton mass compatible with WMAP. On the 
other hand, the hierarchy problem of particle physics is elegantly solved by 
supersymmetry (see e.g. POj), according to which every known particle comes 
with its superpartner, the sparticle. In supersymmetric models the heavy 
singlet neutrinos have scalar partners with similar masses, the sneutrinos, 
whose properties are ideal for playing the role of the infiaton [H^ I34j . 

Superstring theory includes, apart from the fundamental string, other 
extended objects called p-branes. A special class of p-branes are D(irichlet)p- 
branes, where open strings can end. D-brane physics has motivated the 
brane-world idea, which has attracted a lot of interest over the last years. 
In a brane-world scenario our universe is modeled by a 3-brane embedded 
in a five- dimensional bulk spacetime. In the simpest cases, all the standard 
model fields (open string sector) are confined on the brane, while gravity 
(closed string sector) propagates in the bulk. The brane is a hypersurface 
that splits the five- dimensional manifold into two parts and plays the role 
of a boundary of spacetime. Usually the brane is considered to be infinitely 
thin and the matching conditions can be used to relate the bulk dynamics 
to what we observe on the brane. The model first proposed by Randall and 
Sundrum (RS II) ^^ offers a viable alternative to the standard Kaluza-Klein 
treatment of the extra dimensions and together with various extensions has 
been intensively investigated for its cosmological consequences (see e.g. J2S| 
and for reviews |4Uj). 

In four dimensions, the Einstein tensor is the only second-rank tensor 
that (i) is symmetric, (ii) is divergence free, (iii) it depends only on the 
metric and its first derivatives, and (iv) is linear in second derivatives of the 
metric. However, in D > 4 dimensions more complicated tensors with the 
above properties exist. For example, in five dimensions the second order 
Lovelock tensor reads 

TJ D D O D DC o jycd D I Dcde p 

n-ab — J^J^ab — -^-tiac-ttb ~ ^-^ riacbd + J^a J^bcde 

- ^gab{R^ - ^RcdR"" + R'^'^Rcdes) (5.1) 

and can be obtained from an action containing the GB term |1T] 
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^GB = R — ^RabR^ + R°' '^ Rabcd (5.2) 

Higher order curvature terms appear also in the low-energy effective field 
equations arising in string theory. Brane-worlds are string-inspired and so it 
is natural to include such terms in the five-dimensional field equations. 

It is important to note that in the context of extra dimensions and the 
brane-world idea one obtains on the brane a generalized Friedmann equa- 
tion, which is different from the usual one of conventional four-dimensional 
cosmology. This means that the rate of expansion of the universe in this 
novel cosmology is altered and accordingly the description of the physics in 
the early universe can be different from the standard treatment. So it would 
be very interesting to study the cosmological implications of these new ideas 
about extra dimensions and braneworlds. The Friedmann-like equation for 
a GB brane-world has been derived in |^ |331 1^ ■ 

Sneutrino inflation in the context of Randall-Sundrum type II model 
has been analyzed in 03]. However, it would be interesting to study the 
effect of the GB term. After all, this term is a high energy modification to 
general relativity and as such it is expected to be important in the early 
universe. Furthermore, as it has been shown in ^^, the quadratic potential 
V ^ (j)"^ for the infiaton is observationally more favoured when the GB term 
is present. The purpose of the present work is to discuss sneutrino inflation 
in the context of a GB braneworld. 

The present chapter is organized as follows. There are five sections of 
which this introduction is the first. In section 2 we describe sneutrino in- 
flation in a GB brane-world. Section 3 contains the discussion of reheating, 
gravitino production and baryogenesis through leptogenesis. Our results are 
summarized in section 4 and we conclude with a discussion section 5. 



5.2 Sneutrino inflation in a GB brane-world 

5.2.1 GB brane-world 

Here we review GB brane-world, following essentially J3S] • The five- dimensio- 
nal bulk action for the GB braneworld scenario is given by 



S = ^ f (fx^/^^g[-2A, + R 

ZK^ J 

+ a{R'- ARabR"' + RabcdR"'"")] 

(fx^/^\ + Sm^t (5.3) 

brane 
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where a > is the GB couphng, which has dimensions of length'^, A > 
is the brane tension, A5 < is the bulk cosmological constant and S'mat 
denotes the matter action. The fundamental energy scale of gravity is the 
five-dimensional scale M5 with k^ = 87r/M|. For the discussion to follow we 
define a new mass scale through the relation a = 1/M^. 

The GB term may be viewed as the lowest-order stringy correction to 
the five-dimensional Einstein-Hilbert action with a <^ 1/Ai^, where 1/fi is 
the bulk curvature scale, \R\ ~ /i^. The Randall- Sundrum type models are 
recovered for a = 0. Moreover, for an anti-de Sitter bulk, it follows that 



A. 



-3;u2(2-0, where 



^ = Aafi^ < 1 



(5.4) 



Imposing a Z2 reflection symmetry across the brane in an anti-de Sitter 
bulk and assuming that a perfect fluid matter source is confined on the brane, 
one obtains the modified Friedmann equation 



4{p + \) = 2fiJl + 






3 - e + 2e- 






(5.5) 



This can be rewritten in the useful form 

,2 r 



H 



2 ^ A!_ 



:i - cosh 



2X 



1 



(5.6) 



where x is a dimensionless measure of the energy density p on the brane 
defined by 

p -|- A = m^ sinh X (5.7) 



with 



rrir 



8/i^(l- 0- 



311/8 



(5.^ 



the characteristic GB energy scale. 

The requirement that one should recover general relativity at low energies 
leads to the relation 



K^ 



/i 



1+e 



K,^ 



(5.9) 



where k^ = Sn/Mpi and Mpi is the four-dimensional Planck scale. Since 
.^ <^ 1, we have /i ~ M^/M^i. Furthermore, the brane tension is fine-tuned 
to zero effective cosmological constant on the brane 



kI\ = 2/i(3 - 



(5.10) 
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The GB energy scale rua is larger than the RS energy scale A^''^, since we 
consider that the GB term is a correction to RS gravity. Using ()5.1U|) this 
implies jlE] ^ ^ 0.15, which is consistent with Eq. ()5.4|1 . 

Expanding Eq. ()5.6p in x, we find three regimes for the dynamical history 
of the brane universe 



p>m^ 


^ H^^ 


2 2 
LI Kf 


■2 -i 

(GB) 


(5.11) 


m^ > p > A 




(RS) 


(5.12) 


p< A 


=> H' 


^ ^ 
~ 3 p 


(GR) 


(5.13) 



Eqs. (I5.1H) - (J5.13|1 are considerably simpler than the full Friedmann equation 
and for infiation we shall assume the first one (GB). 

5.2.2 Chaotic inflation in a GB brane- world 

We will consider the case in which the energy momentum on the brane is 
dominated by the sneutrino infiaton field confined on the brane with a self- 
interaction potential V{(f)) = (1/2)M^ 0^, where M is the mass of the sneu- 
trino field. The field is a function of time only, as dictated by the isotropy 
and homogeneity of the observed four- dimensional universe. A homogeneous 
scalar field behaves like a perfect fiuid with pressure p = (1/2)0^ — V and 
energy density p = (1/2)0^ + V. We shall assume that there is no energy ex- 
change between the brane and the bulk, so the energy-momentum tensor T^,^ 
of the scalar field is conserved, that is V^T^^ = 0. In terms of the pressure 
p and the energy density p the continuity equation takes the form 

p + 3H{p + p) = (5.14) 

where H is the Hubble parameter H = a/ a. This is equivalent to the equation 
of motion for the scalar field 

+ 3/70 + \/'(0) =0 (5.15) 

the Klein-Gordon equation for in a Robertson- Walker background. The 
equation that governs the dynamics of the expansion of the universe is the 
Friedmann-like equation of the previous subsection. Infiation takes place in 
the early stages of the evolution of the universe, so we suppose that infiation 
takes place in the GB high energy regime 

(2 2 \ 2/3 
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In the slow-roll approximation the slope and the curvature of the potential 
must satisfy the two constraints e <^ 1 and \ri\ <^ 1, where e and r] are the 
two slow-roll parameters which are defined by 

^--^ (5-17) 

V" 
,^^_ (5.18) 

In this approximation the equation of motion for the scalar field takes the 

form 

V 

(5.19) 



3H 
while the generalized Friedmann equation becomes {V ^ 0^) 

H' - (^ Vj (5.20) 

The number of e-folds during inflation is given by 

A^ = In^ = r Hdt (5.21) 



Before presenting all the formulae, it would perhaps be useful at this point 
to describe what follows. Any model of inflation should i) solve the flatness 
and horizon problems, ii) reproduce the amplitude for density perturbations 
(COBE normalization), iii) predict a nearly scale-invariant spectrum, and iv) 
predict very small tensor perturbations. For a strong enough inflation we take 
A^ = 70, which is enough to solve the horizon and flatness problems. Using 
the equations of motion we shall compute the spectral index, as well as the 
scalar and tensor perturbations. We will then flx the remaining parameters 
by requiring that the amplitude of scalar perturbations is reproduced. This 
will lead to a prediction of the spectral index and the tensor-to-scalar ratio. 

According to a recent analysis |17j, at 1 — cr 

A, ~ 2 X 10"^ (5.22) 



-0.048 <ns-l< 0.016 (5.23) 

with As the amplitude of the density perturbations and Us the spectral index. 
On large cosmological scales, data jlTj give for the tensor perturbations 

r < 0.47 95% c.l. (5.24) 
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with r the tensor-to-scalar ratio defined as r = 16 Af/A^ (consistent with the 
normalization of Ref. [23| in the low energy limit), where At is the amplitude 
of the tensor perturbations. 

In the slow-roll approximation the number of e-folds and the slow-roll 
parameters are given by the formulae 



' ^ ^ (5.25) 



fi'^KJ 



9^5/3 

V" f 4e 

3^2/3 



/i2fi;2 



2/3 

(5.26) 



/i-Ztg \ '/' /"^-^ V^2/3 



,2,^2 



where 0end is the value of the infiaton at the end of infiation, which is deter- 
mined from the condition that the maximum of e, \ri\ equals unity, and the 
* denotes the point at which observable quantities are computed. The main 
cosmological constraint (normalization condition) comes from the amplitude 
of the scalar perturbations 



4 H"^ 

A, = ^ -— (5.28) 

where the right-hand side is evaluated at the horizon-crossing when the co- 
moving scale equals the Hubble radius during infiation and Mpi = 1.22 x 
10^^ GeV is the four-dimensional Planck mass. In the present context the 
amplitude of the scalar perturbations is given by 

144\/8/3 / ^2^2X^/3 

The spectral index for the scalar perturbations n^ is given in terms of the 
slow-roll parameters 

d InA^^ 
ns-l = -rn^ = 2r]-6e 5.30 

dink 

and is found to be 

2 AT — 3 

ris = = 0.98 (5.31) 

2N ^ ' 

while the tensor-to-scalar ratio r is given by 

3M2M2 

* "' (5.32) 



2iV3/2jvfM| 
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with At the amphtude of the tensor perturbations 

A = F^ T^- (5.33) 

Sa/tt Mpi 

where again the right-hand side is evaluated at the horizon-crossing. Taking 
the normalization condition into account we obtain for M 

M = 3.4 X 10"^ — ^^ (5.34) 

Ms ^ ^ 



and for the tensor-to-scalar ratio 






r = 75.3 "^^/' (5.35) 



5.3 Reheating, gravitino production and lep- 
togenesis 

5.3.1 Reheating 

We start by introducing three heavy right-handed neutrinos Ni which only 
interact with leptons and Higgs. The superpotential that describes their 
interactions is jlH] 

W = haN.L^H^ (5.36) 

where fia is the matrix for the Yukawa couplings, H^ is the superfield of 
the Higgs doublet that couples to up-type quarks and La (a = e, /i, r) is the 
superfield of the lepton doublets. We assume that the scalar partner of the 
lightest right-handed neutrino plays the role of the infiaton. After inflation 
the infiaton decays into normal particles which quickly thermalize. This is 
the way the universe reenters the radiation dominated era. The sneutrino 
infiaton decays into leptons and Higgs and their antiparticles according to 
the superpotential ()5.36j) and the decay rate is given by [IHI 

r^ = i-/2M (5.37) 

with M the sneutrino mass and p = ^^ |/iaP- The reheating temperature 
after inflation is defined by assuming instantaneous conversion of the infiaton 
energy into radiation, when the decay rate of the infiaton F^ equals the 
expansion rate H. In GB braneworld cosmology H is given by 

H.{^) pV3 (,38, 
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and in the radiation dominated era the energy density of the universe is given 
by 

P = PR = Qeff — T^ (5.39) 

with Qeff = 228.75 the effective number of relativistic degrees of freedom in 
the MSSM for T > 1 TeV. Thus we obtain 



The condition H{Tr) = F,^ gives for the reheating temperature 

After inflation, the direct out-of-equiUbrium decays of the sneutrino inflaton 
generate the lepton asymmetry which is partially converted into a baryon 
asymmetry via sphaleron effects. This requires that Tr < M or that 

I < { 15M| J <"^) 

5.3.2 Gravitino production 

Any viable inflationary model should avoid the gravitino problem jSO] • This 
means that for unstable gravitinos that decay after Big-Bang Nucleosynthesis 
(BBN), their decay products should not alter the abundances of the light 
elements in the universe that BBN predicts. This requirement sets an upper 
bound for the gravitino abundance 

V3/2 = < (5.43) 



n. 



7 



m3/2 



with m3/2 ~ 100 GeV — 1 TeV the gravitino mass, n^ the photon number 
density and (max a parameter related to the maximum gravitino abundance 
allowed by the BBN predictions. According to the analysis of the authors 
of EH, (max = 5 X 10-12 Q^y fo^ ^^^^ ^ ^QO GcV . To find the gravitino 
abundance one has to integrate Boltzmann equation 

^ + 3Hns/2 = C3/2iT) (5.44) 

with 6*3/2(7") the collision term responsible for the thermal production of 
gravitinos as a function of the temperature T < Tr. The rate for the thermal 
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production of gravitinos is dominated by QCD processes since the strong 
coupling is considerably larger than the electroweak couplings. Taking into 
account 10 two-body processes involving left-handed quarks, squarks, gluons 
and gluinos, the authors of [52^ computed the collision term C^/2{T) in the 
framework of supersymmetric QCD. They obtained 

QMT)=a(T)|l+Mr)^j^ (5,45) 

where rrig ~ 1 TeV is the gluino mass and a{T), 6(T) are two slowly- varying 
functions of the temperature, estimated to be jUj 

a{TR) = 2.38, b{Tn) = 0.13 (5.46) 

If we assume that the quantity sa^ is constant during the expansion of 
the universe, where a is the scale factor and s is the entropy density s = 
heff (27r^T'^)/45, then the integration of Boltzmann equation gives 

heff{TR) H[TR)n^[TR) 

with hf.ff the effective number of relativistic degrees of freedom. For T ^ 
1 TeV all particles are relativistic and for the MSSM h^ffiTR) ~ g^jfiTR) = 
915/4 = 228.75, while /ie//(T) = 43/11 for T < 1 MeV. Thus, using flOH|l 
with m3/2 = 100 GeV one is led to the following upper bound for the reheat- 
ing temperature 

Tr < 1.63 X 10-^ "' „,. = To 5.48 

At this point we should also check whether the contribution of the gravitinos 
to the energy density of the universe is compatible with the observed matter 
density of the universe, n^h^ < 0.143 HI], where h = (H/lOO) ^^^^. From 
the gravitino abundance we can calculate their normalized density 

^3/2/1^ = "^3/2 ^3/2^70 ^V«.^ (5-49) 

with n^Q = 3.15 X 10^'^^ GeV^ the photon density today and pcr = 8.07 x 
10~^'^h'^ GeV^ the critical density. For m3/2 = 100 GeV we obtain 



,5/3 
ny^h^ = 1.86 X 10^ ''''^^,,, (5.50) 
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Using the WMAP bound on the matter density of the universe, flmh'^ < 0.143 
we get the following relation between Tr, M5 and M^ 

Tr < 8.54 X 10"^ ^ ^,, — 5.51 

which is less stringent than the constraint ()5.48p coming from BBN. 



5.3.3 Direct leptogenesis from sneutrino decay 

Any lepton asymmetry Y^ = ul/s produced before the electroweak phase 
transition is partially converted into a baryon asymmetry Yb = ris/s via 
sphaleron effects |^. The resulting Yg is 

Yb = CYl (5.52) 

with the fraction C computed to be C = —8/15 in the MSSM |H1]. The lepton 
asymmetry, in turn, is generated by the direct out-of-equilibrium decays of 
the sneutrino infiaton after inflation and is given by 



Vi = ||. (5,53) 

with e the CP asymmetry in the sneutrino decays. For convenience we 
parametrize the CP asymmetry in the form 

e = e""^^ sinSi (5.54) 

where 6^ is an effective leptogenesis phase and e"^"'^ is the maximum asym- 
metry which is given by [55j 

,".» . i. £^M^ (5.55) 

Sir v^stn^p 

with V = 17AGeV the electroweak scale, tanjS the ratio of the vevs of the two 
Higgs doublets of the MSSM and Am^^^ = 2.6 x 10"^ eV"^ the mass squared 
difference measured in atmospheric neutrino oscillation experiments. For 
simplicity we shall take sinf3 ~ 1 (large tan(3 regime), in which case the 
maximum CP asymmetry is given by 

'- = 2x10-° (i54f) (5.56) 
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Combining the above formulae we obtain 



-iii™„x I / ^R 



y^ = 8xlO-^^|.zn<5,|(^^^) (5.57) 



From the WMAP data [TT] we know that 

r/B = — = 6.1 X 10"^° (5.58) 

If we recall that the entropy density for relativistic degrees of freedom is 
s = hf.ff^T^ and that the number density for photons is n^ = -^^T^, one 
easily obtains for today that s = 7.04n^. Thus, using ()5.57|) we have 

T^^l^^GeV (5.59) 

from which we get a lower bound for the rehetaing temperature 

Tr > 1.08 X 10^ GeV (5.60) 



5.4 Results 

Let us summarize the results obtained above. We take M5 and M^, to be two 
independent mass scales, in principle anywhere between the four- dimensional 
Planck mass Mpi and the electroweak scale, v ~ 200 GeV. First we present 
all the constraints that have to be satisfied. We have mentioned that ^ <^ 1 
and that in the GB regime p ^ m\. These lead to the constraints 



2M| 

72" 

and 



M. » ^ (5.61) 



M,<^M (5.62) 

respectively. On the other hand, the sneutrino drives inflation and simulta- 
neously produces the lepton asymmetry through its direct out-of-equilibrium 
decay after the inflationary era. This requires the reheating temperature to 
be smaller than the sneutrino inflaton mass, namely Tr < M. Furthermore, 
the gravitino abundance constraint requires Tr <Tq. So we see that the re- 
heating temperature has to be lower than both M and Tq. Now the question 
arises, whether M is larger than Tq or vise versa. We have checked that for 
M5 and M=K in their allowed range, M is always larger than Tq. Thus, the 

78 



5.4. Results 



requirement that Tr < Tq also guarantees that Tr < M. Hence, for given M5 
and M^,, the reheating temperature is bounded both from below and from 
above as follows 

1.08 X 10^ GeV <Tr<To (5.63) 

Of course, Tq should not be lower than the minimum of the reheating tem- 
perature 

To > 1.08 X 10*^ GeV (5.64) 

Combining all the constraints mentioned above we find an upper bound for 

M, < 3 X 10" GeV (5.65) 

Then, for a given value for M^, M5 has to range between a maximum and a 
minimum value. If M5 gets too small, the tensor-to-scalar ratio gets larger 
than the observed value, while if M5 gets too large, then the constraint ()5.61|) 
or ()5.62|) is not satisfied. For example, for the extremum values of M^ 

• For M, = 3 X 10^1 GeV 

1.31 X 10^^ GeV < M5 < 2.42 x 10^^ GeV (5.66) 

• while for M, = 200 GeV 

9.97 X 10* GeV < M5 < 5.3 x 10^^ GeV (5.67) 

We see that M5 can be very close to the unification scale Mqut ~ 10^^ GeV 
(but remains lower than that) and not lower than 10* GeV . Interestingly, our 
findings are compatible with experiments to probe deviations from Newton's 
law, which currently imply that M5 > 10* GeV \^\. Finally, for all the 
allowed values of M5 and M^,, we find that the constraint ()5.64|) is always 
satisfied and that the tensor perturbations are always negligible. 

So far we have treated M^, as a phenomenological parameter of the model. 
However, the GB coupling a is related to the string mass scale Mstring and 
it is defined to be a = I / {^M'^^^^^g) jM]. Thus, M^ = SM^^^^. M-theory 
seems to allow arbitrary values for the string scale. Experimental limits 
imply that is is not lower than 0{TeV). If the string scale is around a 
few TeV |57| EHl ISH] , observation of novel effects in forthcoming experiments 
becomes a realistic possibility (see e.g. [00] )• For the special case Mgiring = 
7 TeV or M* = 19.81 TeV we obtain 

2.13 X 10^° GeV < Mg < 2.45 x 10^^ GeV (5.68) 
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For the minimuin value M5 = 2.13 x 10^° GeV we obtain for M, tensor-to- 
scalar ratio and reheating temperature the following 

r = marginal (5.69) 

M = 3.28 X 10^2 GeV (5.70) 

and 

1.08 X 10^ GeV <Tr< 4.34 x 10^° GeV (5.71) 

while for the maximum value M5 = 2.45 x 10^'^ GeV we obtain 

r = 3.56 X 10~^ (5.72) 

M = 2.85 X 10^° GeV (5.73) 

and 

1.08 X 10^ GeV <Tn< 6.33 x 10® GeV (5.74) 

Finally, for the Yukawa coupling /^ we find 

• for M5 = 2.13 X 10^° GeV, f < 6.79 x IO-2, 

• while for Mg = 2.45 x lO^^ Q^y ^ p < 553 ^ iQ-e 

However, phenomenological issues such as neutrino masses and axion scale, 
seem more natural if Mgtring is in the range of 10^*^ — 10^^ GeV [HI] centered 
around lO^^ Q^y Yot the case Mstring ~ 10^^ GeV or M, = 3 x lO^^ GeV 
as mentioned already we obtain 

1.31 X 10^^ GeV < Ms < 1.42 x 10^^ GeV (5.75) 

For the minimum value M5 = 1.31 x 10^^ GeV we obtain for M, tensor-to- 
scalar ratio and reheating temperature the following 

r = marginal (5.76) 

M = 3.27 X 10^^ GeV (5.77) 

and 

1.08 X 10^ GeV < Tr < 4.33 x 10^° GeV (5.78) 

while for the maximum value M5 = 1.42 x 10^^ GeV we obtain 

r = 0.4 (5.79) 

M = 3.01 X 10^^ GeV (5.80) 

and 

1.08 X 10^ GeV < Tr < 4.12 x 10^° GeV (5.81) 

Finally, for the Yukawa coupling /^ we find 
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• for Ms = 1.31 X 10^5 GeV, p < 0.07, 

• while for M5 = 1.42 x 10^^ GeV, p < 0.06 

Note that in contrast to the standard four-dimensional [3^ or to the Randall- 
Sundrum sneutrino inflation |13] scenarios, in all cases treated above, the 
Yukawa couphng /^ in the presence of the GB term need not be unnaturally 
small. 



5.5 Conclusions 

In the present work we have examined sneutrino inflation in the GB brane- 
world. The GB term appears in the low-energy effective fleld equations 
of string theories and it is the lowest order stringy correction to the flve- 
dimensional Einstein gravity. Inflation is driven by the sneutrino inflaton, 
which is the scalar superpartner of the lightest of the heavy singlet neutri- 
nos, that might explain in a natural way the tiny neutrino masses via the 
seesaw mechanism. The sneutrino inflaton, apart from driving inflation, also 
produces the lepton asymmetry that partially is converted to the baryon 
asymmetry via sphaleron effects. We flnd that we can get a viable inflation- 
ary model that reproduces the correct amphtude for density perturbations 
and predicts a nearly scale-invariant spectrum and negligible tensor pertur- 
bations. Furthermore, the reheating temperature after inflation is such that 
the gravitino does not upset the BBN results and the required lepton asym- 
metry is generated. Our analysis shows that all these are simultaneously 
achieved for a wide range of values of the flve-dimensional Planck mass M5 
and the mass scale M^, set by the GB coupling. 
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We consider hybrid inflation in the braneworld scenario. In particular, we 
consider inflation in global supersymmetry with the D-terms in the scalar 
potential for the inflaton field to be the dominant ones (D-term inflation). 
We find that D-term dominated inflation can naturally accomodate all re- 
quirements of the successful hybrid inflationary model also in the framework 
of D-brane cosmology with global supersymmetry. The reheating tempera- 
ture after inflation can be high enough (~ 10^° GeV or higher) for successful 
thermal leptogenesis. 



6.1 Introduction 

Recently there has been considerable interest in higher dimensional cosmo- 
logical models. In those models our four-dimensional world lives on a three- 
dimensional extended object (brane) which is embedded in a higher dimen- 
sional space (bulk). The models of this kind are string-inspired ones, as it 
is known that in Type I string theory 162;: there are two sectors, the open 
and the closed ones, and that the theory contains extended objects, called 
D-branes, where open strings can end. The fields in the closed sector (in- 
cluding gravity) can propagate in the bulk, whereas the fields in the open 
sector are confined to the brane. In such string-inspired scenarios the extra 
dimensions need not be small jHZl EB] and in fact they can even be non- 
compact |22] • It is important to note that in the context of extra dimensions 
and the braneworld idea one discovers a generalized Friedmann equation, 
which is different from the usual Friedmann equation in conventional cos- 
mology. This means that the rate of expansion of the universe in this novel 
cosmology is altered and accordingly the physics in the early universe can be 
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different from wliat we know already. So it would be very interesting to study 
tlie cosmological implications of tfiese new ideas about extra dimensions and 
braneworlds. Perliaps tlie best laboratory for such a study is inflation [T5] . 
which has become the standard paradigm in the Big-Bang cosmology and 
which is in favour after the recent discovery from WMAP satellite (see e.g 
|llj ) that the universe is almost fiat. It is known that there is not a the- 
ory for inflation yet. All we have is a big collection of inflationary models. 
The single-field models for inflation, such as 'new' |05| or 'chaotic' [HH], are 
characterized by the disadvantage that they require 'tiny' coupling constants 
in order to reproduce the observational data. This difficulty was overcome 
by Linde who proposed, in the context of non-supersymmetric GUTS, the 
hybrid inflationary scenario jHE]- We remark that before that, the authors in 
|U7j worked out a string-inspired version of hybrid inflation. It turns out that 
one can consider hybrid inflation in supersymmetric theories (for a review on 
supersymmetry and supergravity see fO]) too. In fact, inflation looks more 
natural in supersymmetric theories rather in non-supersymmetric ones [US] . 
In a supersymmetric theory, the tree-level potential is the sum of an F-term 
and a D-term. These two terms have rather different properties and in all 
inflationary models only one of them dominates [On]- The case of F-term 
inflation (where F-terms dominate) was considered for the first time in [THj . 
while the case of D-term inflation (where D-terms dominate) was considered 
in [7T]. In fact, if one considers supergravity then D-term inflation looks more 
promising, since it avoids the problem associated with the inflaton mass [71] . 
F-term inflation in braneworld was studied in [^2]. In the present note we 
discuss the implications of D-term inflation. Before proceeding our discus- 
sion, let us specify our setup. The braneworld model that we shall consider 
is the supersymmetric version of the RS II model (see e.g. [2S])- However, 
the cosmological solution of this extended model is the same as that in the 
non-supersymmetric model, since Einstein's equations belong to the bosonic 
part. The only sourse in the bulk is a five- dimensional cosmological constant. 
There is matter confined to the brane and during inflation, which is the cos- 
mological era we shall be interested in, this matter is dominated by a scalar 
field, called the inflaton field 0. 

The present chapter consists of six sections of which this introduction 
is the first. We present D-term inflation in the second section and brane 
cosmology in the third. Our results for the inflationary dynamics on the 
brane are discussed in the fourth section. We discuss reheating after inflation 
in the fifth section and finally we conclude in the sixth section. 
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6.2 D-term inflation 

In this section we explain what D-term inflation is, following essentially 
Inflation, by definition, breaks global supersymmetry since it requires a non- 
zero cosmological constant V (false vacuum energy of the inflaton). For a 
D-term spontaneous breaking of supersymmetry a term linear in the auxiliary 
field D is needed (Fayet-Iliopoulos mechanism I73|). If the theory contains 
an abelian U{1) gauge symmetry (anomalous or not), the Fayet-Iliopoulos 
D-term 

^[d'eV = ^D (6.1) 



where V is the vector superfield, is supersymmetric and gauge invariant and 
therefore allowed by the symmetries. We remark that an anomalous U{1) 
symmetry is usually present in string theories and the anomaly is cancelled 
by the Green-Schwarz mechanism. However, here we will consider a non- 
anomalous U{1) gauge symmetry. In the context of global supersymmetry, 
D-term inflation is derived from the superpotential 

W = A$<I)+$_ (6.2) 

where $, $_, $+ are three chiral superfields and A is the superpotential cou- 
pling. Under the U{1) gauge symmetry the three chiral superfields have 
sharges Q$ = 0,(5$+ = +1 and Q$_ = —1, respectively. The superpotential 
given above leads to the following expression for the scalar potential 

v(0+,0_, 101) = A2(i0p(i0+p + i0-n + i0+0-n + ^(i0+p - 10-p + 0' 

(6.3) 
where is the scalar component of the superfield $, (f)± are the scalar com- 
ponents of the superfields $-1-, (yf is the gauge coupling of the ^(1) symmetry 
and ^ is a Fayet-Iliopoulos term, chosen to be positive. The global minimum 
is supersymmetry conserving, but the gauge group f/(l) is spontaneously 
broken 

(0) = (0+)=O, (0-) = v^ (6.4) 

However, if we minimize the potential, for fixed values of 0, with respect to 
other fields, we find that for > 0^ = f a/^, the minimum is at 0+ = 0_ = 0. 
Thus, for > 0c and 0+ = 0_ = the tree-level potential has a vanishing 
curvature in the direction and large positive curvature in the remaining 
two directions m^ = A^|0p ± g^C,- 

For arbitrary large the tree-level value of the potential remains constant 
and equal to Vq = {g^/2)^'^, thus plays naturally the role of an inflaton 
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field. Along the inflationary trajectory the F-term vanishes and the universe 
is dominated by the D-term, which splits the masses in the <l>_|_ and $_ 
superfields, resulting to the one-loop effective potential for the infiaton field. 
The radiative corrections are given by the Coleman- Weinberg formula [H] 

^V^-ioo, = ^Y.^-ir^<^n'^ (6-5) 

i 

where A stands for a renormalization scale which does not affect physical 
quantities and the sum extends over all helicity states i, with fermion num- 
ber Fi and mass squared mf. The radiative corrections given by the above 
formula lead to the following effective potential for D-term inflation 

The end of inflation is determined either by the failure of the slow-roll con- 
ditions or when approaches 0c- 



6.3 Effective gravitational equations on the 
brane 

Here we review the basic equations of brane cosmology. We work essentially 
in the context of Randall- Sundrum II model J22]. In the bulk there is just a 
cosmological constant A5, whereas on the brane there is matter with energy- 
momentum tensor r^j,. Also, the brane has a tension T. The five- dimensional 
Planck mass is denoted by M5. If Einstein's equations hold in the five- 
dimensional bulk, then it has been shown in |^ that the effective four- 
dimensional Einstein's equations induced on the brane can be written as 

„ . Svr / 8^ \2 7-1 /^ ^,\ 

G^,u + ^49f^u = ■J^^i^'^ + ij^) ^^^ - Ef,u (6.7) 

p 5 

where g^u is the induced metric on the brane, 7r^,y = u'^ Vi' + 1 3^^'^ '^^P '^"^ ~ 
l^fiaTu ~ ^''"^S'/x!/) ^4 is the effective four-dimensional cosmological constant, 
Mp is the usual four-dimensional Planck mass and E^^, = C^p^ Ua n^ g^ g" is 
a projection of the five-dimensional Weyl tensor Ca/3pa, where n" is the unit 
vector normal to the brane. The tensors vr^i, and Ep_u describe the influence 
of the bulk in brane dynamics. The five- dimensional quantities are related 
to the corresponding four-dimensional ones through the relations 

''- - {IS ^'-'^ 
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and 

in / 47rT2\ ^ ^ 

^' = Mf (^' "" W J "■" 

In a cosmological model in which the induced metric on the brane g^^ has 
the form of a spatially fiat Friedmann-Robertson- Walker model, with scale 
factor a{t), the Friedmann-like equation on the brane has the generalized 
form (see e.g. the second reference in |25j ) 

where C is an integration constant arising from Ef^,^. The cosmological con- 
stant term and the term linear in p are familiar from the four- dimensional 
conventional cosmology. The extra terms, i.e. the "dark radiation" term 
and the term quadratic in p, are there because of the presence of the extra 
dimension. Adopting the Randall-Sundrum fine-tuning 

47rT2 
A» = -3^ (6.11) 

the four-dimensional cosmological constant vanishes. Furthermore, the term 
with the integration constant C will be rapidly diluted during inflation and 
can be ignored. So the generalized Friedmann equation takes the final form 

We notice that in the low density regime p -C T we recover the usual Fried- 
mann equation. However, in the high energy regime p^T the unity can be 
neglected and then the Friedmann-like equation becomes 

Note that in this regime the Hubble parameter is linear in p, while in con- 
ventional cosmology it goes with the square root of p. 



6.4 Inflationary dynamics on the brane 

As already mentioned, we will consider the case in which the energy-mome- 
ntum on the brane is dominated by a scalar field confined on the brane 
with a self- interaction potential V{(j)) given in (jfi.fij) . The field is a function 

86 



6.4. Inflationary dynamics on the brane 



of time only, as dictated by the isotropy and homogeneity of the observed 
four-dimensional universe. A homogeneous scalar field behaves like a perfect 
fluid with pressure p = (1/2)0^ — V and energy density p = (1/2)0^ + V. 
There is no energy exchange between the brane and the bulk, so the energy- 
momentum tensor T^j, of the scalar field is conserved, that is VT^j^^ = 0. 
This is equivalent to the continuity equation for the pressure p and the energy 
density p 

p + 3H{p + p) = (6.14) 

where H is the Hubble parameter H = a/ a. Therefore we get the equation 
of motion for the scalar field 0, which is the following 

+ 3i70 + \/'(0) =0 (6.15) 

This is of course the Klein-Gordon equation for a scalar field in a Robertson- 
Walker background. The equation that governs the dynamics of the expan- 
sion of the universe is the Friedmann-like equation of the previous section. 
Inflation takes place in the early stages of the evolution of the universe, so in 
the Friedmann equation the extra term dominates and therefore the equation 
for the scale factor is 

In the slow-roll approximation the slope and the curvature of the potential 
must satisfy the two constraints e <^ 1 and |?7| ^ 1, where e and t] are the 
two slow-roll parameters which are defined by 

^ - -§-2 (6-17) 

V" 

In this approximation the equation of motion for the scalar field takes the 

form 

V 



^ 3H 
while the generalized Friedmann equation becomes {V ^ 0^) 



(6.19) 



ij2 ^ i^ (g 20) 

3TM2 ^ ^ 



The number of e-folds during inflation is given by 



ttf 






/ 



N = In-^ = I Hdt (6.21) 
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For a strong enough inflation we take A^ = 60. In the slow-roll approximation 
the number of e-folds and the slow-roll parameters are given by the formulae 

m 



^.M^ 


fV'Y AT 


16tt 


\V) V 


-,^« 


(V"\ 2T 


Svr 


[vj V 


Svr 


r<t>f y y 


M2 


4 ^' 2T 



(6.22) 
(6.23) 

N^-^ / ir^dft) (6.24) 

The main cosmological constraint comes from the amplitude of the scalar 
perturbations which is given in this new context by jHI] 

^' ~ 75M6 y^ \2t) ^^-^^^ 

where the right-hand side is evaluated at the horizon-crossing when the co- 
moving scale equals the Hubble radius during inflation. Finally, the spectral 
index for the scalar perturbations is given in terms of the slow-roll parameters 

d InAl 

and the tensor-to-scalar ratio is given by 

A^ T , , 

In what follows we will assume that g ~ 0.5 and that inflation ends at 
0c = {9/^) V^- To make sure that the slow-roll conditions are satisfied we 
impose the constraint 

Also, we have assumed that the potential V is much larger than the brane 
tension T. Therefore another constraint to be satisfied is 



AT 



> 1 (6.29) 



Now that we have written all the necessary formulae, we can proceed to the 
presentation of our results. For arbitrary A it is not possible to satisfy both 
the datum from COBE that Ag = 2 x 10~^ and the slow-roll conditions. For 
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this to happen the superpotential couphng A has to be smaller or equal to 
0.0245 (approximately). Then, for a given value for A, the brane tension 
cannot become arbitrarily large because in that case the constraint that the 
potential should be much larger than the brane tension is not satisfied. We 
find the following upper bound for the brane tension T 

T < 10^^ - 10^^ GeV^ (6.30) 

Now that we have set upper bounds for T and A so that our constraints and 
the data from COBE are satisfied, we can compute the spectral index n^ and 
the tensor-to-scalar ratio r. We find for the spectral index Ug = 0.983 — 0.998 
and for the ratio r ~ 10~^ or lower. 

A detailed analysis shows that for a particular value for A (below the 
upper bound of course) the spectral index does not depend on T and is always 
very close to 1. As A becomes smaller and smaller the spectral index slightly 
increases and gets even closer to 1. Also, in all cases the tensor perturbations 
are negligible. Finally, we find that for the maximum value for the brane 
tension, M^"""'"'^ ~ lO^^Gel/and (v^)^„^ = (3.99-6.74) x lO^^Gel/, whereas 
y/^ becomes smaller as T decreases. We note that according to our analysis A 
a priori can take arbitrarily small values. However, this would be unnatural 
and for that reason we do not consider values for A smaller than 5 x 10~^. In 
that case we find that the values of the infiaton remain safely below Planck 
mass and therefore global supersymmetry is a good approximation. Our 
result differs from a similar study [7S|, in which the discussion leads to the 
conclusion that supergravity corrections are important. 



6.5 Reheating 

Finally, let us turn to the discussion of reheating after inflation and to the 
computation of the reheating temperature Tj^. After slow-roll the infiaton 
decays with a decay rate F and the decay products quickly thermalize. This 
is the way the universe reenters the radiation era of standard Big-Bang cos- 
mology. The reheating temperature Tr is related to two more cosmological 
topics, namely the gravitino problem j^U] and the baryogenesis through lep- 
togenesis. In gravity mediated SUSY breaking models and for an interesting 
range of the gravitino mass, ^713/2 ~ 0.1 — ITeV, if the gravitino is unstable 
it has a long lifetime and decays after the BBN. The decay products destroy 
light elements produced by the BBN and hence the primordial abundance 
of the gravitino is constrained from above to keep the success of the BBN. 
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This leads to an upper bound on the reheating temperature Tr after infla- 
tion, since the abundance of the gravitino is proportional to Tr. A detailed 
analysis derived a stringent upper bound Tr ^ 10^ — 10"^ GeV when gravitino 
has hadronic modes [SH EO! • On the other hand, primordial lepton asymme- 
try is converted to baryon asymmetry [ZZ| in the early universe through the 
"sphaleron" effects of the electroweak gauge theory jSSj . This baryogenesis 
through leptogenesis requires a lower bound on the reheating temperature. 
Leptogenesis can be thermal or non-thermal. For a thermal leptogenesis 
Tij ^ 2 X 10^ GeV [ZH] , whereas for non-thermal leptogenesis T^ ^ 10^ GeV 
[THj . It seems that it is impossible to satisfy both constaints for the reheating 
temperature coming from leptogenesis and the gravitino problem. However, 
the authors of |H21 have showed that in the brane world scenario, that we 
discuss here, it is possible to solve the gravitino problem allowing for the 
reheating temperature to be as high as 10^^ GeV. According to reference jS21 
the gravitino abundance is proportional not to the reheating temperature, as 
is the case in conventional cosmology, but to a transition temperature Tt be- 
tween high temperatures (Tr) and low ones (today's temperature Tq). That 
way the requirement for not over-production of gravitino leads to an upper 
bound for this transition temperature and not for the reheating temperature, 
which can be as high as a satisfactory leptogenesis requires. 

The reheating temperature is given by the formula 

^H= V ^^ 6.31 

where Qeff is the effective number of degrees of freedom at the reheating 
temperature and for the MSSM is geff = ^. Assuming that the infiaton 
decays to the lighest of the three heavy right handed neutrinos ip 

<p -^ i) + i) (6.32) 



the decay rate of the infiaton is 



Svr Vv^ 



(6.33) 



where minfi is the infiaton mass. Mi is the smallest of the three neutrino 
mass eigenvalues and iriinfi > 2Mi. The mass of the infiaton is given in 
terms of the coupling constant g and the Fayet-Iliopoulos parameter ^ by 

rriinfi = V2 g \/^ (6.34) 
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Figure 6.1: Reheating temperature Tr versus the superpotential couphng A 
for Ml = 10^° GeV and v^ = 10^^ GeV . 



If the value of the mass of the hghtest right handed neutrino is Mi = 
10^'^ GeV , which is a representative value, then the reheating temperature 
Tr can be one to two orders of magnitude larger than the right handed neu- 
trino mass, depending on the values of the superpotential coupling A and 
the Fayet-Iliopoulos term ^ (see Figures IFTTl and IFT^ . So we see that the 
reheating temperature is of the right order of magnitude for thermal lepto- 
genesis. When the right handed neutrino mass increases (remaining though 
smaller than minfi/2), the reheating temperature increases too and in fact 
it goes like ~ \/M^. For example, if A = 0.01 and y^ = 10^"^ GeV, then 
Tr = 4.15 X 10^1 GeV for Mi = 6 x 10^ GeV and Tr = 3.79 x 10^^ Q^y fo^. 
Ml = 5 X 10^^ GeV. Finally, for a given Mi and a given ^, when A increases 
then Tr decreases, but only slighty so as to remain of the right order of mag- 
nitude for a successful leptogenesis (see Figure inH}. Also, for a given value 
of Ml and A, when C, increases, Tr increases also (see Figure | 



6.6 Conclusions 

To summarize, we have reexamined supersymmetric D-term dominated hy- 
brid inflation in brane cosmology. We have found that we can reproduce the 
observational data provided that each of the brane tension, five- dimensional 
Planck mass and the superpotential coupling does not exceed a particular 
value. For a given value for the superpotential coupling, when the brane 
tension takes the maximum allowed value then the scale of inflation ^/^ is of 
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Figure 6.2: Reheating temperature Tr versus the Fayet-Ihopoulos term C, for 
superpotential couphng A = 0.01 and Mi = 10^° GeV. 

the order of ~ 10^^ GeV. This value of the inflationary scale is lower than 
the (supersymmetric) GUT scale, but close to it. Also, we have found that 
for natural values of the superpotential coupling A the inflaton field cannot 
take large values and stays well below the four-dimensional Planck mass, 
consistent with the global supersymmetry approximation adopted here. Fur- 
thermore, we have seen that our results are compatible with the correspond- 
ing results in the standard four- dimensional cosmology. This means that the 
advantages of the hybrid model are naturally preserved in the framework of 
brane cosmology. Finally, our study shows that the reheating temperature 
after inflation can naturally be of order lO^'' GeV (or larger) allowing for a 
successful thermal leptogenesis. 
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Chapter 7 

Brane-Bulk energy exchange 



The role of brane-bulk energy exchange and of an induced gravity term on 
a single braneworld of negative tension and vanishing effective cosmological 
constant is studied. It is shown that for the physically interesting cases of 
dust and radiation a unique global attractor which can realize our present 
universe (accelerating and < Qmo < 1) exists for a wide range of the param- 
eters of the model. For Q^no = 0.3, independently of the other parameters, 
the model predicts that the equation of state for the dark energy today is 
u!de,o = —1-4, while Qmo = 0.03 leads to wde,o = —1-03. In addition, during 
its evolution, wde crosses the wde = — 1 line to smaller values. 



7.1 Introduction 

In cosmologies where the present universe is realized as a finite point during 
the cosmic evolution, the answer to the coincidence question "why it is that 
today Qmo and QoEfi are of the same order of magnitude" , relies on appro- 
priate choice of initial conditions. By contrast, in a scenario in which the 
present universe is in its asymptotic era (close to a fixed point) the answer to 
the above question reduces to an appropriate choice of the parameters of the 
model. However, this latter situation is not easily realized if today's universe 
is accelerating, because: 

If the energy density of a perfect fiuid with equation of state w > —1/3 
of any cosmological system is conserved, all fixed points of the system with 
Qm 7^ are decelerating. 

Indeed, with p the energy density of the perfect fiuid with conservation 
equation p + 3{l + w)Hp = 0, the Hubble equation of an arbitrary cosmology 

93 



Brane-Bulk energy exchange 



can be written in the form 

H^ = 2-i{p + pde) (7.1) 

where 7 = AtiGn/^- Then, the equation governing poE can always be 
brought into the form p^E + 3(1 + wde)Hpde = 0, where wde is time- 
dependent and distiguishes one model from the other. It can be easily seen 
that diytm/^DE) / d\\ia = 3{Qrn/^DE){'WDE — w) and 2q = 1 + 3{wflrn + 
wde^de), where Qm = '^IP/H'^, ^de = ^'JPde/H'^ and q = —d/aH"^. At 
the fixed point (denoted by *) d{Qrn/^DE)/d\na = 0. For Qm* 7^ one 
obtains wde* = w, and 2g=K = 1 + 3w7 > 0. 

Thus, independently of the cosmological model, the only way our accel- 
erating universe with Q^* 7^ can be close to a late time fixed point is by 
violating the standard conservation equation of matter. In 4-dimensional the- 
ories, an accelerating late time cosmological phase characterized by a frozen 
ratio of dark matter/dark energy appears in coupled dark energy scenar- 
ios [HD] as a result of the interaction of the dark matter with other energy- 
momentum components, such as scalar fields. In higher dimensional theories, 
where the universe is represented as a 3-brane, this violation could be the 
result of energy exchange between the brane and the bulk. In particular 
in five dimensions, a universe with fixed points characterized by Qm* 7^ 0, 
g* < was realized in |23 in the context of the Randall-Sundrum braneworld 
scenario with energy influx from the bulk. However, these fixed points can- 
not represent the present universe, since they have Qm* > 2. In this paper 
we present a brane-bulk energy exchange model with induced gravity whose 
global attractor can represent today's universe. 

Let us consider an arbitrary cosmology in the form (jT.ip . Instances of 
such cosmologies arise in braneworld models or in theories with modified 4- 
dimensional actions leading to H^ = f{p), or in cosmologies where poE is 
due to additional fields. Assuming that as a result of some interaction p is 
not conserved, it will satisfy an equation of the form 

p + 3{l + w)Hp = -T (7.2) 

Then, the equation governing poE can always be brought into the form 

Pde + 3(1 + wde)Hpde = T (7.3) 

where wde is time and model dependent. Whenever a fixed point of the 
system satisfies 

H,T, ^0 , p = pDE = (7.4) 
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one obtains 

1 + w 
wde* = -1 - 1 _ (7.5) 

Equation ()7.5j) is model-independent, in the sense that it does not depend 
on the form of T or the function wnEit). For fi^* < 1 equation ()7.5|) gives 
wde* < —1- Specifically, ior w = and Qm* = ^cdm = 0.3 one obtains 
wde* = -1-4, while for fi^* = ^bar = 0.03, wde* = -1-03. 

The cosmology discussed in the present paper has a global attractor of the 
form ()7.4|1 . ()7.5j) [81 . Moreover, the universe during its evolution crosses the 
Wde = — 1 barrier from higher values. This behavior is favored by several 
recent model-independent [82j as well as model-dependent [HHl IHH lEHl tHEj 
analyses of the astronomical data. 



7.2 The model 

We consider the model described by the gravitational brane-bulk action jHZ| 

S= f <fx^^ {M^R - A) + /" d^xV^ {m^R - V) (7.6) 

where i?, R are the Ricci scalars of the bulk metric qab and the induced metric 
hAB = Qab — nAUB respectively {n^ is the unit vector normal to the brane 
and A,B = 0, 1, 2, 3, 5). The bulk cosmological constant is A/2M^ < 0, the 
brane tension is V, and the induced-gravity crossover scale is fc = m? /M^. 

We assume the cosmological bulk ansatz 

ds"^ = -n{t, yfdt^ + a{t, y^-fijdx'dx^ + b{t, yfdy"^ (7.7) 

where 7jj is a maximally symmetric 3-dimensional metric, parametrized by 
the spatial curvature k = —1,0,1. The non-zero components of the five- 
dimensional Einstein tensor are 



a /a b\ n 



.a\a bJ h^ 



a" a! (o! b' 



a a \ a b 



kn 



a2 



(7.8) 



a^ (a' /a' 2n'\ b' fn' 2a' \ 2a" n" 
o^laVa n J b \n a J a n 
a (2h a\ 2a b (h 2a\ b' 



a ( a / Zn a\ Za o /n Za\ o"i , , , 

+ -7^,- )-- + I )-l\- ^^^^ ^-9 

77,^ la\n aJ a b\n a) b) 

(7.10) 
(7.11) 
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where primes indicate derivatives with respect to y, while dots derivatives 

with respect to t. The five- dimensional Einstein equations take the usual 

form 

1 



G 



AC 



2M3 



T, 



AC\tot 



where 



rpA\ 

-l-Cltot 



T'A I I rpA I 

C \v,B 'T~ -'- C \m,B 



T'A I _[_ n^A\ _[_ n^A\ 



is the total energy-momentum tensor, 



(7.12) 
(7.13) 



rpA\ 

J-C \v,B 



rpA\ 

-'r 

J-C l"i,b 



-C \v,( 
nA\ 



diag(-A, -A, -A, -A, -A) 
diag(-y,-r,-r,-V,0)^ 



diag(-p,p,p,p,0) 



Ky) 



(7.14) 
(7.15) 

(7.16) 



T^\m.B is any possible additional energy- momentum in the bulk, the brane 
matter content T^\m,b consists of a perfect fiuid with energy density p and 
pressure p, while the contributions arising from the scalar curvature of the 
brane are given by 



J- n 



\ind 



J j I ind 



6m /d kn\6{y) 



-^ + 



2m (a 



■w 



;,2 



a 
2ah 

an 



^2d^ fcn^\^,%) 
a a? ) ^ h 



(7.17) 
(7.18) 



Assuming a Z2 symmetry around the brane, the singular part of equations 
()7.12|1 gives the matching conditions 



p + F 



n' 



Unhr, 



ttoho 

2p + 3p-V 
12M3 



12M3 ^ 2^2 Va2 ^ 



d^ 



kn'i 



+ 



2nl 



2dn 



al 



2aoho 



knt 



(7.19) 
(7.20) 



(the subscript o denotes the value on the brane), while from the 05, 55 
components of equations ()7.12j) we obtain 



a' 



tin ^ tXn ilr. 



a' 



nf^Oo ^ afoK _ ^ 

iIq (Iq Gjq Oq Gjq 



Tr 



05 



Tli Ittr. 



Qo 



+ 



an /a. 



an \a, 



rir. 



Ur. 



khl 



n. - Abl 



6M3 



(7.21) 
(7.22) 
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where Tq^,T5^ are the 05 and 55 components of TAc\m,B evaluated on the 
brane. Substituting the expressions ()7.19|) . ()7.20p in equations ()7.2H) . ()7.22)] . 
we obtain the semi- conservation law and the Raychaudhuri equation 












r^^{p + 3p-2V)i r^^{p + 3p-2V){p + V) 



+ 






24m2 

,2 



+ 



1 



^te + 4 



+ 



144m4 

rlip + V) 
12m2 






(7.23) 



(7.24) 



where Ho = ao/aoUo is the Hubble parameter of the brane. One can easily 
check that in the limit m — > 0, equation ()7.24|) reduces to the corresponding 
second order equation of the model without R |23| . Energy exchange between 
the brane and the bulk has also been investigated in jHHl lEHl EOj • 

Since only the 55 component oiTAc\m,B enters equation ()7.24|1 . one can 
derive a cosmological system that is largely independent of the bulk dynam- 
ics, if at the position of the brane the contribution of this component relative 
to the bulk vacuum energy is much less important than the brane matter 
relative to the brane vacuum energy, or schematically 



T^5 

-^5 



A 



< 



(7.25) 



Then, for |A| not much larger than the Randall-Sundrum value V"^ / 12M^ , 
the term T^ in equation ()7.24|) can be ignored. Alternatively, the term T^ 
can be ignored in equation ()7.24|1 if simply 






A 



<1 



(7.26) 



Note that relations ()7.25|1 and ()7.2(ij) are only boundary conditions for T^, 
which in a realistic description in terms of bulk fields will be translated 
into boundary conditions on these fields. In the special case where fl7.25p . 
()7.26|) are valid throughout the bulk, the latter remains unperturbed by the 
exchange of energy with the brane. 

One can now check that a first integral of equation ()7.24|1 is 



Ht 



2H'„ 



p + V 6 

\ 



3k\ 

+ 



p 



V k 



Qim? 
4 / A 



+ 



r2Vl2M3 



X 
3tI 



(7.27) 
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with X satisfying 



X + ^rioHoX 



rjnjT 
m?ho 



Hi 



p + V k 
\ 



(7.28) 



and T = 2T^ is the discontinuity across the brane of the 05 component of 
the bulk energy-momentum tensor. The solution of ()7.27|) for Ho is 



and equation ()7.28|) becomes 

2nl T f r, 



2(P + V) 12 



r/i^ 



A 
M3 



+ X 



(7.29) 



m?ho 



2v^L 



2(P + V) 12 



771^ 



A 



+ X 



(7.30) 



At this point we find it convenient to employ a coordinate frame in which 
ho = rio = 1 va. the above equations. This can be achieved by using Gauss 
normal coordinates with 6(t, z) = 1, and by going to the temporal gauge on 
the brane with Uo = 1. It is also convenient to define the parameters 



A 

/^ 

7 



12 
2 



2V_ 

V 

Qim? 

1 

12m2 
1 



A 

Jp 



(7.31) 

(7.32) 

(7.33) 
(7.34) 



For a perfect fluid on the brane with equation of state p = wp our system is 
described by equations ()7.23|1 . ()7.29j) . ()7.3()|1 . which simplify to (we omit the 
subscript o in the following) 



H' 



p + 3il + w)Hp = -T 

}^ 

p + 2^p± I3y/X + 247P + X - -ij 



X 



+ AHx = 247t(i ± ^VA + 247p + 



X 



while the second order equation ()7.24|) for the scale factor becomes 

A + 6(l-3w)7p 



a 
a 



p-{l + 3w)^p±/3- 



VA + 247P + X 



(7.35) 
(7.36) 

(7.37) 
(7.38) 
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Finally, setting ip = ^A + 24-fp + x, equations dZ^SI), dZSZI), (TT^ take the 
form 

H^ = fi + 2-fp±(3tlj-— (7.39) 

,l. + 2H{,!.-^±^^l^^^)^±?2f (7.40) 

°=,-(l + 3u.hp± / + "' "'"'"" (7.41) 

Throughout, we will assume T{p) = Ap", with u > 0, A constant parameters 
j2Zl 1211 • Notice that the system of equations ()7.35|) - ()7.37|) has the influx- 
outflow symmetry T -^ —T, H — > — -ff, t — » —t. For T = the system 
reduces to the cosmology studied in j9l] . 

We will be referring to the upper (lower) ± solution as Branch A (Branch 
B). We shall be interested in a model that reduces to the Randall-Sundrum 
vacuum in the absence of matter, i.e. it has vanishing effective cosmological 
constant. This is achieved for p = =p/3vA, which, given that m^V + 12M^ 
is negative (positive) for branches A (B), is equivalent to the fine-tuning 
A = —V'^/12M^. Notice that for Branch A, V is necessarily negative. Cos- 
mologies with negative brane tension in the induced gravity scenario have 
also been discussed in [021 • 

Consider the case k = 0. The system possesses the obvious fixed point 
{p^,H^,il)^) = (0,0, a/A). However, for sgn{H)T < there are non-trivial 
fixed points, which are found by setting p = ^ = in equations ()7.35|] . 
(fTini) . For w <l/3 these are: 

^/?^ \9 9 = 2/i + (1 - 3w)7p, 

9{1 + wypl ^ ^ "^ (7.42) 

± v/9(l + w)272p* + 4/92[A -I- 6(1 - 3w)-fp,] 

H. = - ^.^[^*\ (7.43) 

3(1 + w)p^ 

iPl ± ^^i^TP^V-* - [A + 6(1 - 3w)7P*] = (7.44) 

Equation ()7.4H) gives 



which is positive, and also, it has the same form (as a function of p^) as in 
the absence of R. The deceleration parameter is found to have the value 

g. = -1 (7.46) 
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which means H^, = 0. Furthermore, at this fixed point we find 



a„. . HI^ . M±i^.pr- (7.47) 



Equation ()7.42|) . when expressed in terms of flm*, has only one root for each 
branch 



/3 6(1 - 3w)/3 ± V\{1 -3w- 4f]-i~ 
2^ (2fi^l + l + 3M;)(f]^l-l) 



/>* = T^ ' 7oo 1 V o\7n i~T; (7-48) 



However, it can be seen from ()7.48p that for — 1 < u' < 1/3 and Q^n* < 1 the 
Branch B is inconsistent with equation ()7.42|) . On the contrary. Branch A 
with —l<w< 1/3 and ^m* < 1 is consistent for < 6(1 — 3w)l3 + vA(l — 
3i(; — 4^2"^) < 3a/4(1 — 3w)^/3^ ~ (1 + w^X. Thus, since we are interested in 
reahzing the present universe as a fixed point. Branch B should be rejected, 
and from now on we will only consider Branch A. So, we have seen until 
now that for negative brane tension, we can have a fixed point of our model 
with acceleration and < Qm* < 1- This behavior is qualitatively different 
from the one obtained in the context of the model presented in |2ZI (for 
— 1/3 < w < 1/3), where for positive brane tension we have Qm* > 2, while 
for negative brane tension the universe necessarily exhibited deceleration; 
therefore, in that model the idea that the present universe is close to a fixed 
point could not be realized. 

Concerning the negative brane tension the following remarks are in or- 
der: (a) In the conventional, non-supersymmetric setting, it is well known 
that a negative tension brane with or without induced gravity is accompa- 
nied by tachyonic bulk gravitational modes ^22j; however, including the GB 
corrections relevant at high-energies, the tachyonic modes can be completely 
removed for a suitable range of the parameters [HI], (b) As shown in ^^, 
in supersymmetric theories, spacetimes with two branes of opposite tension 
are stable; in particular, there is no instability due to expanding "balooning" 
modes on the negative brane. It is, however, unclear what happens in models 
with supersymmetry unbroken in the bulk but softly broken on the brane. 
(c) Finally, it has been shown jHEl that with appropriate choice of boundary 
conditions, both at the linearized level as well as in the full theory, the grav- 
itational potential of a mass on a negative tension brane has the correct 1/r 
attractive behaviour. 
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7.3 Critical point analysis 

We shall restrict ourselves to the flat case k = 0. In order to study the 
dynamics of the system, it is convenient to use (dimensionless) flatness pa- 
rameters such that the state space is compact [HZj. Defining 



Ur. 



27P 

D2 






H 
D 



where D = sjll'^ — /x, we obtain the equations 



A /|/i|c(;„ 



;i + 3w)(u;m- 1)^- 
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2Z(1 



Z' 



^^ 1 - Z^ - 3{1 - 3w)/3^fi-^iUr, 

1 — U^m 
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1 — iOm. 



'l-Z')i' 



(7.49) 

(7.50) 
(7.51) 



l + 3w 



-UJr. 



(7.52) 



with ' = d/dr = D ^d/dt. Note that — 1 < ^ < 1, while both tu's satisfy 
< a; < 1. The deceleration parameter is given by 



Q 



1 
Z2 



l + 3w 



.OJr, 



-UJ„ 



3(1 - 3w)/?V^^« 



1 



iOr. 



(7.53) 



and H' = —HZ{q +1). The system of equations ()7.5H) - ()7.52j) inherits from 
equations ()7.35|l - ()7.37|l the symmetry A — *• —A, Z -^ —Z, r -^ — r. The 
system written in the new variables contains only three parameters. However, 
going back to the physical quantities H, p one will need specific values of 
two more parameters. 

It is obvious that the points with \Z\ = 1 have H = oo. Therefore, 
from ()7.39|) it arises that the infinite density p = oo big bang (big crunch) 
singularity, when it appears, is represented by one of the points with Z = 1 
[Z = —1). The points with Um = 1, \Z\ 7^ 1,0 have u'^ = 00, Z' = 00 and 
finite p, H; for w < 1/3, one has in addition d/a = +00, i.e. divergent 4D 
curvature scalar on the brane. 

The system possesses, generically, the fixed point (a) {u!m*,^tp*, Z^,) = 
(0,1,0), which corresponds to the fixed point {p^, H^jijj^,) = (0,0, V^) dis- 
cussed above. For u < 3/2 there are in addition the fixed points (b) (a;m*, o;^*, 
Z*) = (0,1,1) and (c) {uJm*,^-4!*, Z^:) = (0,1,-1). All these critical points 
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z/<3/2 


z/ = 3/2 


z/>3/2 


No. of F.P. 


1 


or 1 


1 


Nature 


A 


A 


S 



Table 7.1: The fixed points for w = 0, influx. 

are either non-hyperbolic, or their characteristic matrix is not defined at 
all; thus, their stability cannot be studied by first order perturbation anal- 
ysis. In cases like these, one can find non-conventional behaviors (such as 
saddle-nodes and cusps jUH]) of the flow-chart near the critical points. There 



are two more candidate fixed points (d) 



iUn 



, ^Ip^ 



1,0,1) and (e) 



(a;m*,i^V*' ^*) ~ (1)0,-1), whose existence cannot be confirmed directly 
from the dynamical system, since they make equations ()7.51|) . ()7.52|) unde- 
termined. Apart from the above, there are other critical points given by 



A / |/i|^-n 



v-l 



m + w)z. 



'\-Z^- 



1-^(1 



UJn 



2[l - {I -Zt 



2\2l 



(7.54) 



(7.55) 



They exist only for AZ^ < and correspond to the ones given by equa- 
tions ()7.42j) - ()7.44|l . For the physically interesting case w = with influx we 
scanned the parameter space and were convinced that for u ^ 3/2 there is 
always only one fixed point; for u < 3/2 this is an attractor (A), while for 
u > 3/2 this is a saddle (S). For w = 0, u = 3/2 there is either one fixed point 
(attractor) or no fixed points, depending on the parameters. For the other 
characteristic value w = 1/3, we concluded that for u < 3/2 there is only 
one fixed point (attractor), for z/ > 2 there is only one fixed point (saddle), 
while for 3/2 < u < 2 there are either two fixed points (one attractor and one 
saddle) or no fixed points at all, depending on the parameters. For w = 1/3, 
u = 3/2 there is either one fixed point (attractor) or no fixed points. Finally, 
for w = 1/3, z/ = 2 there is either one fixed point (saddle) or no fixed points. 
These results were obtained numerically for a wide range of parameters and 
are summarized in Tables mi and 17^ 

The approach to an attractor described by the linear approximation of 
()7.51|) - ()7.52|) is exponential in r and takes infinite time r for the universe to 
reach it. Given that near this fixed point the relation between the cosmic 
time t and the time r is linear, we conclude that it also takes infinite cosmic 
time to reach the attractor. 
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z/<3/2 


z/ = 3/2 


3/2 < z/ < 2 


z/ = 2 


z/ > 2 


No. of F.P. 


1 


or 1 


Oor 2 


or 1 


1 


Nature 


A 


A 


A,S 


s 


S 



Table 7.2: The fixed points for w = 1/3, influx. 

Defining e = sgn{H), we see from fl7.5HI - ()7.52|l that the hues Z = e 
{u < 3/2), uJm = are orbits of the system. Furthermore, the family of solu- 
tions with Z ^ e and dZ/duJm = Z' /u'^ «i is approximately described for 
u < 3/2 by uj'^ = e(l + 3w)uJm{^m — 1), and thus, they move away from the 



point {uJm*,Z^] 



%l) 



while they approach the point 



.UJr, 



Zjfi 



= (1,-1)- 

In addition, the solution of this equation is uJrn = [l + ce''*^^+^"'^'^] ^, with c > 
an integration constant. Using this solution in equation H'/H = —Z{q + 1) 
we find that for w = 1/3, H/ Ho = ^/uJ^/{l — ujm), where Ho is another 
integration constant. Then, the equation for uJm{t) becomes dujm/dt = 
—2euJm\/ Hoidm ^ /^(l ~ "-^m)^, and can be integrated giving t as a function 
of uJm or H. Therefore, in the region of the big bang/big crunch singularity 
one obtains a(t) ~ Vet, p{t) ~ t~^, as in the standard radiation dominated 
big-bang scenario. This means that for v < 3/2 the energy exchange has no 
observable effects close to the big bang/big crunch singularity. 

Since our proposal relies on the existence of an attractor, we shall restrict 
ourselves to the case z/ < 3/2. It is convenient to discuss the four possible 
cases separately: 

(i) w = with influx. The generic behavior of the solutions of equations 
()7.5H) - ()7.52|) is shown in Figure 17.11 We see that all the expanding 
solutions approach the global attractor. Furthermore, there is a class 
of collapsing solutions which bounce to expanding ones. Finally, there 
are solutions which collapse all during their lifetime to a state with 
finite p and H. The physically interesting solutions are those in the 
upper part of the diagram emanating from the big bang {uj, Z) ^ (1,1)- 
These solutions start with a period of deceleration. The subsequent 
evolution depends on the value of 3/5^/|yu|, which determines the relative 
position of the dashed and dotted lines. Specifically, for 3/5^/|/i| > 1 
(the case of FigureEH]) one distinguishes two possible classes of universe 
evolution. In the first, the universe crosses the dashed line entering the 
acceleration era still with wde > — 1, and finally it crosses the dotted 
line to wde < — 1 approaching the attractor. In the second, while in the 
deceleration era, it first crosses the dotted line to wde < — 1, and then 
the dashed line entering the eternally accelerating era. For 3/5^/|yu| < 1, 
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Figure 7.1: Influx, w = 0, u < 3/2. The arrows show the direction of 
increasing cosmic time. The dotted hne corresponds to wde = —1. The 
region inside (outside) the dashed hne corresponds to acceleration (deceler- 
ation). The region with Z > represents expansion, while Z < represents 
collapse. The present universe is supposed to be close to the global attractor. 

the dotted line lies above the dashed line, and, consequently, only the 
second class of trajectories exists. To connect with the discussion in the 
introduction, notice that the Friedmann equation ()7.39|1 can be written 
in the form ()7.H) with dark energy p^E = {Pip + /^)/27- Using ()7.4U|) . 
the equation for pde takes the form ()7.3p with 



-1 



wde 



3(1 — Ur, 



2Z'^ -LJm-l- 6(1 - 3w) 



jj, Z'^ -Um 



(7.56) 



The global attractor ()7.42|1 - ()7.44|1 satisfies relations ()7.4j) and conse- 
quently, wde evolves to the value wde* given by ()7.5|1 . As for the 
bouncing solutions, they approach the attractor after they cross the 



line Z 



UJ„ 



where wde jumps from +oo to — oo; however, the evo- 



lution of the observable quantities is regular. 

(ii) w = with outflow. The generic behavior in this case is obtained from 
Figure mi bv the substitution Z — > —Z and r — > — r, which reflects the 
diagram with respect to the uj^ axis and converts attractors to repelers. 

(iii) w = 1/3 with outflow. Figure 17^ depicts the flow diagram of this case. 
Even though in the case of radiation in general wde > —1/3 from 
equation ()7.56|1 . there are both acceleration and deceleration regions. 
Furthermore, from equation ()7.5|1 it is Qm* > 1- 
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z 1 




Figure 7.2: Outflow, w = 1/3, u < 3/2. The arrows show the direction 
of increasing cosmic time. The region inside (outside) the dashed line cor- 
responds to acceleration (deceleration). The region with Z > represents 
expansion, while Z < represents collapse. 

(iv) w = 1/3 with influx. This arises like in (ii) by reflection of Figure 17^ 
and resembles Figure mi 



7.4 Conclusions 

In this work, we studied the role of brane-bulk energy exchange on the cosmo- 
logical evolution of a brane with negative tension, zero effective cosmological 
constant, and in the presence of the induced curvature scalar term in the 
action. Adopting the physically motivated p^ power-law form for the energy 
transfer and assuming a cosmological constant in the bulk, an autonomous 
system of equations was isolated. In this scenario, the "dark energy" is a 
result of the geometry and the brane-bulk energy exchange. The negative 
tension of the brane is necessary in order to realize the present universe (ac- 
celerating with < Qrno < 1) as being close to a future fixed point of the 
evolution equations. We studied the possible cosmologies using bounded nor- 
malized variables and the corresponding global phase portraits were obtained. 
By studying the number and nature of the fixed points we demonstrated nu- 
mericaly that our present universe can be easily realized as a late-time fixed 
point of the evolution. This provides an alternative answer to the coinci- 
dence problem in cosmology, which does not require specific fine-tuning of 
the initial data. Furthermore, the equation of state for the dark energy at the 
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attractor is uniquely specified by the value fimo- Remarkably, for fi^o = 0.3, 
one obtains wdea = —1-4, independently of the other parameters, while for 
the other suggestive value Qmo = 0.03, wdea = —1.03. In the past, the 
function wde crosses the line wde = — 1 to larger values. 

It would be interesting to investigate if the above partial success of the 
present scenario persists after one tries to fit the supernova data and the 
detailed CMB spectum |99j. Of course, the nature of the content of the bulk 
and of the mechanism of energy exchange with the brane is another crucial 
open question, which we hope to deal with in a future publication. 
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Conclusions 



We have studied dark matter, inflation and dark energy in the brane-world 
scenario. These are three topics of fundamental importance for modern cos- 
mology. According to recent observations (Cosmic Microwave Background, 
Supernovae, galaxy serveys), most of the energy content of our universe is in 
the form of dark matter and dark energy. In addition, inflation is responsible 
for producing the cosmological fluctuations for the formation of the struc- 
ture that we observe today. The brane-world idea is inspired from M/string 
theory and although brane models are not yet derived from the fundamen- 
tal theory, at least they contain the key features of string theory, such as 
extra dimensions, higher-dimensional objects (branes), higher order correc- 
tions to gravity (Gauss-Bonnet) etc. Inflation and the dark sector have been 
addressed in the framework of standard four-dimensional cosmology. How- 
ever, it is challenging to try to study them using alternative gravitational 
theories such as braneworlds. Furthermore, since string theory claims to 
give us a fundamental description of nature, it is important to study what 
kind of cosmology they predict. The essence of the brane-world scenario 
is that the standard model, with its matter and gauge interactions, is lo- 
calized on a three-dimensional hypersurface, called brane, embedded in a 
higher-dimensional spacetime, called the bulk. Gravitons, the mediators of 
the gravitational interaction, are free to propagate into the whole bulk and 
thus connect the standard model sector with the dynamics of the extra di- 
mensions. Brane models are capable of giving non-conventional cosmologies. 
In standard four-dimensional cosmology the Hubble parameter goes like the 
square of the energy density. However, in brane models this is not true any 
more. The relation between the Hubble parameter and the energy density 
is more complicated and the specific form depends on the model at hand. 
These novel cosmologies can be used in a two-fold way: On one hand they 
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help us in attacking in a different context longstanding problems of standard 
cosmology. On the other hand they offer us an ideal laboratory for testing 
ideas coming from M/string theory. We believe that brane cosmology is an 
exciting field and that it deserves further study. 
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